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1 Introduction 



At TeV colliders, electroweak radiative corrections are strongly enhanced by loga- 
rithms of the type \n{Q'^/M^i) [1^4]. These logarithmic corrections affect every reaction 
that involves electroweakly interacting particles and is characterized by scattering ener- 
gies Q Mw,z- The impact of the enhanced electroweak corrections at high-energy 
colliders has been investigated both in complete one-loop calculations and in logarithmic 
approximation for several specific reactions, including gauge-boson pair production at the 
ILC [4,5] and the LHC [6], gauge-boson scattering [2,7], fermion-pair production in e+e~ 
collisions [3,8], Drell-Yan processes at the LHC [9], heavy-quark production [10], single- 
gauge-boson plus jet production at the LHC [11], Higgs production in vector-boson fusion 
at the LHC [12], and three-jet production in e+e" coUisions [13]. The results of these 
studies demonstrate that at energies Q ^ 1 TeV the size of the electroweak corrections 
can reach, depending on the process, up to tens of per cent at one loop and several per 
cent at two loops. 

At high energies, the dominant effects can be described in a systematic way by treating 
the electroweak corrections in the asymptotic limit Q/Mw,z oo, where all masses of 
order Mw are formally handled as infinitesimally small parameters. In this limit, the 
electroweak corrections appear as a sequence of logarithms of the form a' In"* (Q'^/ 2) > 
with j <2l, which diverge asymptotically. These logarithms have a two-fold origin. The 
renormalization of ultraviolet (UV) singularities at the scale /xr < Mw,z yields terms 
of the form ln((5^///^) and, in addition, the interactions of the initial- and final-state 
particles with soft and/or coUinear gauge bosons give rise to ln((5^/M^ 2) terms that 
represent mass singularities. 

As pointed out in Ref. [14], soft and collinear electroweak logarithms are present not 
only in those physical observables that are exclusive with respect to real radiation of Z 
and W bosons, but even in fully inclusive observables [14, 15] . However, there is no need to 
combine corrections resulting from virtual and real Z/W bosons in the same observable. 
Thus, in this paper we will consider only exclusive observables, which do not include real 
Z/W-boson emission and contain only virtual contributions. 

Since they originate from UV, soft, and collinear singularities, electroweak logarithmic 
corrections have universal properties that can be studied in a process-independent way 
and reveal interesting analogies between QED, QCD, and electroweak interactions. At 
one loop, the leading logarithms (LLs) and next-to-leading logarithms (NLLs) factorize 
and are described by a general formula that applies to arbitrary Standard-Model pro- 
cesses [16, 17]. The properties of electroweak logarithmic corrections beyond one loop can 
be investigated by means of infrared evolution equations that describe the electroweak 
interactions in terms of two regimes corresponding to SU(2)xU(l) and Ucm(l) symmet- 
ric gauge theories [18-22]. This approach makes use of resummation techniques that 
were derived in the context of symmetric gauge theories (QED, QCD), thereby assuming 
that electroweak symmetry-breaking effects other than the splitting into two symmetric 
regimes are negligible in next-to-next-to-leading logarithmic (NNLL) approximation. An 
alternative method, which also relies on factorization and exponentiation of the loga- 
rithmic corrections, is based on soft-collinear effective theory [23]. Using this method, 
ln{Q^/M^2) corrections to the Sudakov form factor for massless and massive fermions 
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have been calculated and used to compute electroweak corrections to four-fermion pro- 
cesses at the LHC [24,25]. Alternatively, the electroweak logarithms can be extracted 
explicitly from Feynman diagrams at two loops [26-32]. Calculations of this type explic- 
itly implement all aspects of electroweak symmetry breaking through the Feynman rules 
and thus provide a strong check of the results obtained via evolution equations or effective 
theories. 

So far, all existing diagrammatic results are in agreement with the resummation pre- 
scriptions. However, up to now only a small subset of logarithms and processes has been 
computed explicitly at two loops: while the LLs [26] and the angular-dependent subset 
of the NLLs [27] have been derived for arbitrary processes, complete diagrammatic cal- 
culations at (or beyond) the NLL level exist only for matrix elements involving massless 
external fermions [28-32]. In the literature, no explicit 2-loop NLL calculation exists for 
reactions involving massive scattering particles. 

In this paper we derive the two-loop NLL corrections for general n-fermion processes 
/1/2 /s • • • /n involving an arbitrary number of massless and massive fermions, i.e. 
leptons, light or heavy quarks. We consider the limit where all kinematical invariants 
are of order » -^w,Z) the top mass — as well as the Higgs mass — is of the same 
order as Mw,z- Apart from the top quark all other fermions, including bottom quarks, 
are treated as massless particles. Soft and coUinear singularities from virtual photons are 
regularized dimensionally and arise as e-poles in D = 4 — 2e dimensions. For consistency, 
the same power counting is applied to ln((5^/Afw) and 1/e singularities. Thus, in NLL 
approximation we include all e"'' ln-'~'^((5^/M^) terms with total power j = 2, 1 at one 
loop and j = 4, 3 at two loops. We explicitly show that the photonic singularities can be 
factorized in a gauge-invariant electromagnetic term, in such a way that the remaining 
part of the corrections — which is finite, gauge invariant, and does not depend on the 
scheme adopted to regularize photonic singularities — contains only ln{Q^/M^) terms. 
The divergences contained in the electromagnetic term cancel if real-photon emission is 
included. 

We utilize the technique that we introduced in Ref. [32] to derive the NLL two-loop 
corrections for massless n-fermion processes. This method is based on collincar Ward 
identities which permit to factorize the soft-coUinear contributions from the n-fermion 
tree-level amplitude and isolate them in process-independent two-loop integrals. The lat- 
ter are evaluated to NLL accuracy using an automatized algorithm based on the sector- 
decomposition technique [33] and, alternatively, the method of expansion by regions com- 
bined with Mellin-Barnes representations (see Ref. [31] and references therein). 

The treatment of n-fermion processes that involve massive top quarks implies two new 
aspects: additional Feynman diagrams resulting from Yukawa interactions, and top-mass 
terms that render the loop integrals more involved. The latter is the one that requires more 
effort from the calculational point of view, since every topology has to be evaluated for 
several different combinations of massive/massless internal and external lines. This part 
of the calculation was automatized in the framework of the above-mentioned algorithms, 
thereby doing an important step towards a complete NLL analysis of all processes that 
involve massive particles. A brief summary of the results presented here was anticipated 
in Ref. [34]. 
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The paper is organized as follows. Section [2] contains definitions and conventions used 
in the calculation. In Sect. [3] we review the techniques used to extract UV and mass 
singularities for the case of massless fermion scattering [32] and extend them to the case 
of massive fermions. The one-loop counterterms and explicit results for the renormalized 
one- loop amplitude are given in Sect. |H The complete two- loop results are presented 
in Sect. |5] including a discussion of the Yukawa-coupling contributions and the two-loop 
renormalization. Section [6] is devoted to the discussion of our final results. Explicit 
results for the various contributing one- and two-loop diagrams are presented in App. \M 
Further appendices contain the definition of the loop integrals (App. [B]) and relations 
between them (App. [C]). Specific results for four-particle processes involving external 
fermions and gluons, including a comparison with the results of Refs. [24,25], can be 
found in App. [Dl 



2 Definitions and conventions 



The calculation is based on the formalism introduced in Ref. [32] for massless fermion 
scattering. Here we summarize the most important conventions and introduce new defi- 
nitions that are needed to describe massive fermions. For more details we refer to Sect. 2 
of Ref. [32]. 

We consider a generic n — process involving an even number n of polarized fermionic 
particles, 

ipi{pi) . . . iPniPn) 0. (2.1) 

The symbols ipi represent n/2 antifermions and n/2 fermions: ipi = for i = 1, . . . ,n/2 
and = for i = n/2 + 1, . . . ,n. The indices = R, L and cxj characterize the chirality 
(see below) and the fermion type (/o-, = Uey^, . . . ,T,u,d, . . . ,t), respectively. All external 
momenta are incoming and on shell, pi = ml. Apart from the top quarks, all other 
fermions (including bottom quarks) are treated as massless particles. 
The matrix element for the process (12.11) reads 



M 



ipi...ifi„ 



n/2 
i=l 



G^^-^^{pu 



! P ri j 



n ^iPj^i^j) 

i=n/2+l 



(2.2) 



where G-i is the corresponding truncated Green function. The spinors fulfil the Dirac 
equation. 



m)u{p, k) = 0, 



-|- m)v{p, k) 



0, 



(2.3) 



and the argument k denotes their chirality. More precisely, the polarization states k = R, L 
are defined in such a way that in the massless limit the spinors are eigenstates of the chiral 
projectors 



1/. 5 

UJR = ujl = -(1 +7 



= = ^(1-7^), 



(2.4) 



I.e. 



ujpu{p, k) = Sf,pU{p, k) + 



m 



ujpv{p, k) = 8^,pv{p, k) + 



m 



(2.5) 
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for m/p^ <^ 1. While massless spinors are exact eigenstates of the chiral projectors, in 
the massive case the spinors are constructed by means of hehcity projectors 

^^R = ^^L = ^(1 + 7V), l^L = ^^R = ^(1 - 7V) (2.6) 



where 



r 



_ (\^ ^0 P 



±- IpI,/^ , for sign(/) = ±l, (2.7) 



m \ ' \p\ 
with (rp) = 0, = —1 and 

fipu(p, k) = 5^pu{p, k), (lpv{p, k) = 5^pv{p, k). (2.8) 

For m/p^ —>■ 0, we have — > p'^/m, and the spinors satisfy (12. 5p in the high-energy hmit. 
Note that in the case of antifermions, chirahty and hehcity play the opposite role, and we 
always use chirality to label the spinors. 

The amplitudes for physical scattering processes, i.e. 2 n — 2 reactions, are easily 
obtained from our results for — > reactions using crossing symmetry. 

2.1 Perturbative and asymptotic expansions 

For the perturbative expansion of matrix elements we write 

M-Im,, M,^(^)'m., (2.9) 

where a = e^/(47r) and in D = 4 — 2e dimensions we include in the definition of a 
normalization factor depending on e, the scale /id of dimensional regularization, and the 
characteristic energy Q of the scattering process. 

The electroweak corrections are evaluated in the region where all kinematical invari- 
ants, rj,,,fc = {pj + - ■ .+Pk)'^, are much larger than the squared masses of the heavy particles 
that enter the loops. 



rj...k\ 



^ M^^ml ^ M^. (2.10) 



In this region, the electroweak corrections are dominated by mass-singular logarithms. 



L = In ( -^1 , (2.11) 

and logarithms of UV origin. Mass singularities that originate from soft and collinear 
massless photons and UV singularities are regularized dimensionally and give rise to 1/e 
poles. The amplitudes are computed as series in L and e, classifying the terms 
according to the total power m of logarithms L and 1/e poles. At / loops, terms with 
m = 2/,2/ — 1, . . . are denoted as leading logarithms (LLs), next-to-leading logarithms 
(NLLs), and so on. The calculation is performed in NLL approximation expanding the 
one- and two-loop terms up to order and e°, respectively. 



-Ml ='E E -Mi,„.,ne"L-+", Y: M,^^^^e-L"^^r (2.12) 

m=l n=— m m=3 n=—m 
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Since the loop corrections depend on various masses, Mw ~ Mz ~ mt ~ Mu, and different 
invariants rjk, the coefficients M.imn (l^-l^P involve logarithms of typcl 

^. = ln(^), ^^•^ = 1^(^^^) fo^^'^^- (2-13) 
For convenience we also define Ijj = and 

L. = >" (I) . (2.14) 

To distinguish terms associated with massless and massive fermions we use the symbols 





rrii = mt 






mi = 1 




mi = 




[0, 






mi = rrit j 



kt = , <^.o = • (2.15) 

\^0, rrii = j \^0, rrii = nit j 

Mass-suppressed corrections of order M^/Q"^ are systematically neglected. 

2.2 Gauge and Yukawa couplings 

The generators associated with the gauge bosons V = A, Z, W"^ are related to the weak 
isospin generators T*, the hypercharge Y, and the electromagnetic charge Q through 

e/^* = ^{t'± iT^) , e/^ = -eQ = -g2S^T^ - gicj^, 
el' = g^c^T' - gts^^ = - e^Q, (2.16) 



where Cw = cos^„ and = sin6'„ denote the sine and cosine of the weak mixing angle 
and gi and g2 are the coupling constants associated with the U(l) and SU(2) groups 
respectively. The generators fl2.16p obey the commutation relations 



jVtjV, 



m E e^'^'^^'u''^ , (2.17) 
V3=A,z,iy± 



where V denotes the complex conjugate of V and the e-tensor is defined in Ref. [32]. The 
SU(2)xU(l) Casimir operator reads 

E ^^^^ = 4f^)' + 4c. with C = E(TT. (2.18) 



V=A,Z,W^ c \ ^ / c -^-^ 

The gauge-boson masses, the vacuum expectation value f , and the couplings fulfil 

^w± = \g2V, Mz = T^g2V, Cw = c^giY^ = s^g2, (2.19) 

where is the hypercharge of the Higgs doublet. Leaving Y$ as a free parameter and 
identifying e = CwS'i, the Gell-Mann-Nishijima relation, which determines the hyper- 
charges of the fermions, reads Q = Y/2 + Y^T^. 



^ This dependence only appears at the next-to- leading level, i.e. for m = 2/ — 1. 
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The Feynman rules for the vector-boson-fermion-antifermion vertices read 




(2.20) 



where IJ^^jk denote the SU(2)xU(l) generators in the fundamental (k = L) or trivial 
[k = R) representation. The chiral projectors Uk in (12.201) can easily be shifted along the 
fermionic lines using anticommutation relationq^ until they meet the spinor of an external 
fermion or antifermion and can be eliminated using (12. 5p . 

It is convenient to adopt a notation that describes the interactions of fermions and 
antifermions in a generic way. To this end, for a generic incoming particle ipi = f^' or 
we define 

I^"^' = \ • (2.21) 

[jy .^ = -iy for = 

With this notation the interaction of a gauge boson V with incoming fermions and an- 
tifermions yields 

:ieYC',,MPi,i^i) and v{pi, Ki)ieYC^,i (2.22) 



Similar expressions are obtained for multiple gauge-boson interactions. Apart from the 
spinors and the reversed order of the Dirac matrices, these expressions differ only in the 
sign of the mass terms in the Dirac propagator^. In practice we perform the calculations 
assuming that certain legs . . .) are incoming fermions and we find that the results 
are independent of the sign of the mass terms in the propagators (only squared mass 
terms give rise to unsuppressed contributions). Thus, the results are directly applicable 
to fermions and antifermions. 

The Feynman rules for the scalar-fermion-antifermion vertices read 



i (wrGjl jR + ^lGJrjl ) , (2.23) 



^ In the high-energy NLL approximation diagrams involving chiral anomalies are irrelevant. Thus we 
can use {7'', 7^} = in _D = 4 — 2e dimensions. 

■^For antifermions the negative sign in the coupling (|2.2ip compensates the negative sign (of slashed 
momenta) due to the opposite fermion and momentum flow. This yields — — i + nii) = (^i + ^ — rrii). 
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where $ = H,x,4>^- Unitarity implies the following relation between left- and right- 
handed coupling matrices: 

We consider only contributions proportional to the top-quark Yukawa coupling, 

and neglect the bottom-quark mass. The non-vanishing components of the right-handed 
coupling matrix read 

Ct^tH = ^tH^ ^ ^bH^ ^ (2.26) 

In analogy with (12.211) . we define 

. G% ,R for ipi = ff 

Gl„, = { , (2.27) 

'G* p=-G* . for ^. = /T 

and similarly for R L. With this notation the interaction of a boson $ with incoming 
fermions = and antifermions (fi = yields similar expressions as in the case of 
gauge interactions f l2.22p . 



In the case of light (ant i) fermions {(pi =lepton or quark of the first two generations), the 
subsequent interaction with gauge bosons Vi, V2, V3 yields coupling factors / 

where the representation of all generators corresponds to the chirality Ki given by the 
spinor of the external particle cpi [see (12.211) ]. For quarks of the third generation (yjj =top 
or bottom), due to Yukawa interactions and top- mass terms in the fermion propagators, 
also couplings with opposite chirality contribute. These are denoted as 



fv _ jV 



(2.29) 



For instance, the subsequent emission of scalars and gauge bosons $i,V2,$3 along a 
heavy-quark line yields coupling factors G*?,, ,,J^^, 

In our results the matrix elements fl2.2|) are often abbreviated as 

M = M'^'-'^\ (2.30) 

and when they are multiplied by gauge- and Yukawa- coupling matrices we write 

MGf'I^'=yM^'-<-^-Glh,I^f^„, etc. (2.31) 



Similar shorthands are used for the couphng matrices fl2.29p . Global gauge invariance 
implies the relation 

ilGt^-Gt^ll= E gViI^^ (2.32) 

between combinations of gauge and Yukawa couplings as well as the charge- conservation 
identity 

n 

MY.ll = 0, (2.33) 

fe=i 

which is fulfilled up to mass-suppressed terms in the high-energy limit. In (12.321) . 
denotes the SU(2) xU(l) generators for the Higgs doublet, which enter the gauge couplings 
of the Higgs boson. 

3 Treatment of ultraviolet and mass singularities 

Large logarithms and 1/e poles originate from UV and mass singularities. These 
contributions can be extracted from one- and two-loop Feynman diagrams within the 
't Hooft-Feynman gauge, using the technique introduced in Ref. [32]. In Sect. 13. ll we 
review this method for the case of massless fermion scattering. The new aspects that 
emerge in the presence of massive fermions are discussed in Sect. 13. 2[ Finally, in Sect. 13.31 
we report on a calculation of the fermionic form factor as a check of the validity of our 
methods. 

3.1 Massless fermions 

Here we give a concise summary of the method presented in Ref. [32]. For a detailed 
discussion we refer to the original paper. 

3.1.1 Origin of mass singularities 

Mass singularities appear in loop diagrams involving soft and / or coUinear gauge bosons 
that couple to external legs. At one loop, the mass singularities of the n-fermion amplitude 
originate from diagrams of the type 

V 

Vl f^ - , (3.1) 

where an electroweak gauge boson, V = A,Z, iy±, couples to one of the external fermions 
(i = 1, ... ,77.) and to either (i) another external fermion or (ii) an internal propagator. 
While diagrams of type (ii) produce only single logarithms of coUinear origin, diagrams 
of type (i) give rise to single and double logarithms. The latter originate from the region 
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where the gauge-boson momentum is soft and colhnear to one of the external fermions. 
At two loops, the following five types of diagrams give rise to NLL mass singularities 






(3.2) 

Here the NLL mass singularities originate from the regions where the gauge boson Vi is 
simultaneously soft and collinear, and the gauge bosons V2 and V3 in the first four types 
of diagrams are soft and/or collinear. 



3.1.2 Factorizable and non-factorizable contributions 



The soft-collinear NLL contributions resulting from the diagrams (13. ip and (13.21) are 
split into factorizable and non-factorizable ones. The one- and two-loop factorizable (F) 
parts result from those diagrams where the virtual gauge bosons couple only to external 
lines. Including sums over gauge bosons and external legs we have 



Ml 



EE 



and 



i=i V=A,Z,W± 




(3.3) 



-M^ = EE E 



i = l j = l Vrr,=A,Z.W± 



+ 




+ E 



fe=i 




1 

6 




k=l 1=1 




(3.4) 

The limit xp^ in (13.31) and (13. 4p indicates that the above diagrams are evaluated 

in the approximation where each of the four-momenta of the various gauge bosons 
is collinear to one of the momenta of the external legs or soft. Where relevant, also 
the contributions of hard regions are taken into account (see Sect. 13. IT^ . The label 
F in the n-fermion tree subdiagrams in (13. 3p and (13.41) indicates that, by definition, the 
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factorizable contributions include only those parts of the above diagrams that are obtained 
by performing the loop integration with the momenta of the gauge bosons set to 
zero in the n-fermion tree subdiagrams. 

These contributions are called factorizable since their logarithmic soft-collinear singu- 
larities factorize from the n-fermion tree amplitude (see below). The remaining contribu- 
tions are called non- factorizable (NF). They comprise those diagrams of type (13. ip and 
(13. 2p that involve gauge bosons Vi coupling to internal (hard) propagators, and the non- 
factorizable parts of the diagrams in (13.31) and (13.41) . In Ref. [32], using coUinear Ward 
identities, it was explicitly shown that all non-factorizable one- and two-loop contributions 
cancel at the amplitude level. 

3.1.3 Soft-collinear approximation 

The soft-collinear singularities are extracted from the above Feynman diagrams using 
the following soft-collinear approximation for the interactions of virtual gauge bosons 
Vi . . .Vn with an incoming (anti)fermion line i 



-2e/^ (p, + QnT- -2e/r {p^ + qiY' 

{Pi + qnY {Pi + qiY 

(3.5) 

Here are the loop momenta of the gauge bosons Vk, and qk = qi + ■ ■ ■ + qk- In practice, 
the coupling of each soft /collinear gauge boson gives rise to a factor —2elY''{pi + quY- 
When applied to the factorizable diagrams in (I3.3l) - (l3.4p . this approximation removes all 
Dirac matrices occuring along the fermionic lines and yields factorized expressions of the 
form 

Ml^ = M^K^,2. (3.6) 

where the NLL corrections -^'1,2, involving coupling factors and logarithmically divergent 
loop integrals, factorize from the n-fermion tree amplitude M.q. 

As discussed in Ref. [32], the soft-collinear approximation (13.51) provides a correct 
description of the gauge-boson-fermion couplings in all soft /collinear regions that are 
relevant for an NLL analysis at one and two loops, with the following exception: The 
soft-collinear approximation is not applicable to topologies of type 




(3.7) 



which contain UV singularities associated with a subdiagram with characteristic scale 
/^foop ^ (^66 next section) or give rise to power sing ularities of 0(1/M2). For these 
diagrams, the soft-collinear approximation (13. 5p can be applied only to the vertices that 
occur outside the one- loop subdiagrams that are depicted as grey blobs in (13. 7p . whereas 
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for the vertices and propagators inside the one-loop sub diagrams we have to apply the 
usual Feynman rules. 

To bring these diagrams in the factorized form fl3.6p . we have to simplify the Dirac 
matrices along the fermionic lines where the soft-coUinear approximation is not applicable. 
To this end we utilize trace projectors that reduce these Dirac matrices to scalar quantities. 
In practice, for the second and third diagram in (13.71) . the Dirac matrices along the line i 
are eliminated using [32] 

Xu{pi, Ki) = J2 ^p^ij i^pX) u{Pi, Ki) = Uij (uj^^X) u{pi, Ki), (3.8) 
p 

where X represents the diagram without jXZ) external spinors, and the projector is 
defined as 

% (r) = Tr [r^,^,] (3.9) 

if rrii = rrij = 0. For a detailed discussion we refer to Sect. 3.2 of Ref. [32]. 

3.1.4 Treatment of ultraviolet singularities 

In addition to logarithmic soft-coUinear singularities, also the logarithms originating 
from UV singularities need to be taken into account. Moreover, the soft-coUinear ap- 
proximation (13.51) can give rise to fake UV logarithms that must carefully be avoided or 
subtracted. 

Since UV singularities produce only a single logarithm per loop, in NLL approximation 
we need to consider only UV-divergent diagrams of one-loop order and their insertion in 
one-loop diagrams (13.31) involving leading soft-collinear singularities. The UV logarithms 
originate from the UV cancellations between bare diagrams and counterterms as 

bare diagrams counterterms 

where /xioop is the characteristic scale of the UV-singular loop diagram, /ir is the renor- 
malization scale, and the term (/^d/Q^)'' — that we always absorb in [see (12. 9p ] — is 
factorized. In our calculation, the UV poles of bare loop diagrams and counterterms are 
removed by means of a minimal subtraction at the scale Q^. As a result, the UV divergent 
terms take the form 

(*n;[(S'-i-H(S)'-ii' 

. ' ' 

bare diagrams counterterms 

which is obviously equivalent to (I3.10p . The advantage of this subtraction is that for 
all bare (sub)diagrams with fifo^^ ~ the UV singularities e~^[((5^//ik)op)'' "~ M do not 
produce large logarithms and are thus negligible. This holds also for fake UV singu- 
larities resulting from the soft-collinear approximation. Thus, apart from the counter- 
terms, we must consider only those (subtracted) UV contributions that originate from 
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bare (sub) diagrams with /i^^op <^ Q^. At one loop this condition is never reahzed in 
the high-energy hmit (12.101) . and in practice we need to consider only the two- loop UV 
contributions which result from the diagrams in (13. 3p through insertion of UV-divergent 
subdiagrams in the lines that are not hard (yufoop ^ Q^)- These contributions corre- 
spond to the diagrams depicted in (13.71) . for which we use — instead of the soft-collinear 
approximation — the projectors (I3.8l) - (l3.9p . thereby ensuring a correct description of the 
UV regions. 

3.2 Massive fermions 

The method outlined in Sect. 13. H is applicable also to processes involving external top 
and/or bottom quarks. However, in the presence of massive fermions, two new aspects 
must be taken into account: mass terms in the top-quark propagators and new diagrams 
resulting from Yukawa interactions. 

3.2.1 Top-mass terms 

Let us recall that, in the high-energy limit (I2.10p . the top-quark mass is treated as 
a small parameter and terms of (9(NLL x m^/Q) are systematically neglected. In this 
approximation, only the following types of mt-terms must be considered: 

(i) Yukawa couplings proportional to mt/iWw; 

(ii) mt-terms acting as coUinear regulators in the denominator of propagators; 

(iii) mt-terms in the numerator of loop integrals of (9(NLL/mt) or (9(NLL/Mw,z,h); i-G- 
integrals that involve power singularities. 

The majority of the mt-terms occuring in one- and two-loop diagrams does not belong 
to these categories and can be set to zero. In particular, by explicit inspection of the 
relevant one- and two-loop integrals, we find that mt-terms of type (iii) occur only inside 
the one-loop insertions in the diagrams (13.70 . Apart from these special cases all other 
mt-terms in the numerators of loop integrals can be set to zero. 

3.2.2 NLL contributions from gauge interactions 

As in the case of massless fermion scattering, the diagrams (I3.ip - (l3.2p with external-leg 
interactions of soft / collinear gauge bosons produce logarithmic mass singularities. These 
diagrams are split into factorizable and non-factorizable parts as discussed in Sect. 13.1.21 
and, in the presence of massive fermions, the Ward identities that are responsible for the 
cancellation of the non-factorizable contributions [32] receive only negligible corrections 
of 0{mt/Q). Thus, as in the case of massless fermions, the non-factorizable parts do not 
contribute and we can restrict ourselves to the calculation of the factorizable parts (13.30 
and M . 

The soft-collinear approximation (13. 5p can easily be adapted to massive fermions by 
including the mt-terms in the denominator of the top-quark propagators and leaving the 
numerator as in the massless case. In principle mt-terms modify also the numerator of 
(13. 5p . However, as observed above, such terms are relevant only in the one- loop insertion 
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diagrams in (13. 71) . Here, as discussed in Sect. 13.1.31 we employ the usual Feynman rules — 
and not the soft-collinear approximation — due to the presence of UV divergences or power 
singularities of C(l/M^). 

To bring the diagrams of type (13.71) in the factorized form (13.61) . we utilize projections 
analogous to (13.81) . In the case of massive fermions, using the two projectors 



and 



we obtain 



n., (r) 



n., (r) 



Tr 



r(^i + mA 



PiPj 



Tr r(^i + mi) 1 - 



PiPj . 



(3.12) 



(3.13) 



(3.14) 



Here, in contrast to the massless case, the Dirac matrices are not projected out completely 
and a term proportional to Ilij{ujpX)'^j remains, which cannot be cast in the factorized 
form (13. 6p . However, after explicit evaluation of the loop integration, we find that these 
Hjj-terms are suppressed in the high-energy limit and only the factorizable terms associ- 
ated with the projector Hj^ contribute. 

The projectors (I3.12p - (l3.14p are applicable in the case where particle i is an incoming 
fermion. For antifermions similar projectors can be constructed, which differ only in the 
sign of the m^-terms. As discussed in Sect. 12. 2[ this difference is irrelevant since only 
squared mass terms produce NLL contributions. Thus, the results derived for fermions 
are directly applicable to antifermions. 

All logarithms of UV origin — also in the presence of Yukawa interactions — are treated 
by means of a minimal subtraction at the scale as explained in Sect. I3.1.4[ 

3.2.3 NLL contributions from Yukawa interactions 

The interaction of external bottom or top legs with scalar bosons is suppressed in all 
soft/collinear regions. Thus, Yukawa interactions can produce NLL contributions only 
through the counterterms (see Sect. 14. 1| ) and the topologies of type (13.70 . which contain 
UV singularities or power singularities of (9(1/M^). In addition to the diagrams of type 
(13. 7p that are present in the massless case, the following new bare diagrams must be taken 
into account: 



EE E E 

i=l J = l Vrr,=A,Z,W± *i=H,X,</>± 



+ 




+ 




(3.15) 
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By explicit inspection we find that the last two diagrams in (13.151) are suppressed in NLL 
accuracy since the VV^ couplings of 0(Mw,z) are not compensated by 1/M terms from 
the loop integrals. The NLL terms resulting from the first three diagrams in fl3.15l) are 
worked out in Sect. 15. 1[ 

3.3 Form- factor checks 

When evaluating the factorizable contributions of arbitrary n-fermion processes, we 
have to eliminate the Dirac matrices along the fermionic lines in order to separate the 
loop integrals from the tree amplitude A^g- For this purpose we use the soft-coUinear 
approximation for gauge interactions presented in Sect. 13. 1.31 and the fact that the Yukawa 
interactions are suppressed in all soft /collinear regions (see Sect. 13.2.31) . 

We have checked the validity of this procedure for the case of the form factor which 
couples a fermion-antifermion pair to an external Abelian field. In our approach, the 
radiative corrections to this form factor are given by the factorizable contributions (13.31) . 
(13. 4p . and (13.151) with n = 2 external fermions (in this case the diagrams with three 
and four legs in (13. 4p do not contribute). Alternatively, we have calculated the one- and 
two-loop corrections to the form factor in the high-energy limit including all Feynman 
diagrams, without neglecting mt-terms in the numerator, and employing projection tech- 
niques (see e.g. Ref. [35]) instead of the soft-collinear approximation. After performing 
the minimal subtraction of the UV divergences as explained in Sect. 13.1.4^ the results of 
the complete form-factor calculation agree with the ones resulting from the factorizable 
contributions (13. 3p . (13. 4p . and (I3.15P in soft-collinear approximation. We have also verified 
by explicit evaluation that, as expected from the discussion in Sect. 13.2.31 the form-factor 
diagrams with Yukawa interactions that are not included in (13.150 are either suppressed or 
vanish in NLL accuracy. Two of these additional Yukawa diagram^ yield non-vanishing 
NLL contributions which are, however, of pure UV origin and are completely removed by 
the minimal UV subtraction. 

4 One-loop results 

The one-loop amplitude gets contributions from the factorizable one-loop diagrams 
(13. 3p and from counterterms, 

M^ = Ml + Mf^. (4.1) 

The results for the loop diagrams are summarized in App. lA.li Here we discuss the one- 
loop renormalization and list the final results for the renormalized one-loop amplitude. 
As shown in Ref. [16] the one-loop NLL corrections factorize, i.e. they can be expressed 
through correction factors that multiply the Born amplitude. Moreover, they can be split 
in a symmetric-electroweak part, an electromagnetic part, which in particular contains 
all soft /collinear singularities associated with photons, and an Mz-dependent part that 
results from the difference in the W- and Z-boson masses. This splitting permits to 
separate the soft /collinear singularities resulting from photons in a gauge-invariant way 
and is very important in view of the discussion of the two-loop contributions. 

"^These are the one-loop diagram (jA.2|) with the gauge boson Vi replaced by a scalar boson and the 
two-loop diagram (jA.7P with the gauge boson V2 replaced by a scalar boson. 
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4.1 One-loop renormalization 



As discussed in Ref . [32] , mass renormalization is not relevant in NLL approximation. 
Thus, counterterm contributions result only from the renormalization of the electroweak 
gauge couplings and the top-quark Yukawa coupling, 



9k,0 



9k + Yl 



1=1 



477 



X (0 

OQk, eo = e 



a 

4:7!- J 



V, 



(4.2) 



and from the renormalization constants associated with the wave functions of the external 
fermions i = 1, . . . ,n, 



1=1 



(4.3) 



All couplings are renormalized in the MS scheme, but we moreover subtract the UV 
singularities both in the bare and the counterterm contributions as explained in Sect. I3.1.4[ 
The counterterm for the Yukawa coupling can be determined from the divergent parts 
of the counterterms to g2, M^^ and using tree-level relations. Assuming that the 
renormalization sea 1^ /iR is of the order of or larger than Mw, we find for the counterterms 



^9, 



(1) NLL Qk 1 , (1) 



5\ 



2 e 

(1) NLL 

2 e 



At 



^1 -1 



^1 -1 



. (1) NLL _e 1,(1) 

2e ^ 



^1 -1 



(4.4) 



where the dependence on the factor {Q"^ / finY is due to the normalization of the expansion 
parameter in (14.21) . and the one- loop /3-function coefficients in the electroweak Standard 
Model (Y^ = 1) are given by 



(1) 



41 
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Qcl' Qsl' "3' ^' Asl ' Ucl 



11 



9 17 

+ 



■ (4.5) 



For later convenience we also define the QED /5-function coefficient, which is determined 
by the light-fermion contributions only. 



^QED 



80 



(4.6) 



where represents the colour factor, i.e. A*"/ = 1 for leptons and = 3 for quarks. 
The renormalization of the mixing parameters Cw and Sw can be determined via (12.191) 
from the renormalization of the coupling constants. 

The contribution of the counterterms (14.41) can easily be absorbed into the Born am- 
plitude 



Mr 



gk=ak{Q^),e=e{Q^),\t=\t{Q^) 



(4.7) 



""We do not identify the renormalization scale /zr and the scale of dimensional regularization /id- 
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2\ NLL 2/ 2 
= e Ut, 





10' -1 




[ An e 




}=eVI){l- 


\ Att e 


i(S)'-: 


} = A?(/i^)|l- 




via the running couplings 

(4.8) 

For practical applications one can use MS input parameters at the scale /iR = Mz, or 
alternatively the on-shell input parameters a(Mz), Mz and Mw, or the input scheme. 
These different schemes are based on input parameters at the electroweak scale and are 
thus equivalent in NLL approximation. In the following we express all one- and two-loop 
results in terms of the Born amplitude at the scale Q^. The notation A^q(Q^) emphasizes 
the fact that the Born amplitude implicitly depends on logarithms ln{Q'^/fi^) via the 
running of the couplings. 

In this setup, the only one-loop NLL counterterm contribution arises from the on- 
shell wave-function renormalization constants dZ^^^ for the fermionic external legs. All 
legs receive contributions from massive weak bosons, whereas the photonic contribution to 
the counterterm for massless external fermions vanishes owing to a cancellation between 
UV and mass singularities within dimensional regularization. The legs with top or bottom 
quarks additionally get Yukawa contributions. After subtraction of the UV poles we find 



(1) NLL 




VtV 



in 



Ml, 



At A 




+ z 



Y 



2e2 



(4.9) 



The Kronecker symbol 5j^t is defined in f l2.15p . Compared with the massless case two 
extra terms appear, one related to the massive top quark in photonic diagrams and one 
originating from the Yukawa couplings. The Yukawa factors are obtained from 



and read 



EA$+^<I) _ Y \ 2 
Lrj (jTj — 



1, for left-handed third-generation quarks = t^,t^,b^,lr, 

2, for right-handed top quarks, ipi = t^,t^, 
0, otherwise. 



(4.10) 



(4.11) 



Finally, the one-loop counterterm for a process with n external fermions in NLL approx- 
imation is obtained as 



M7 = M,iQ') ± \5Z„ 

i=l 



(1) 



(4.12) 



16 



4.2 Renormalized one-loop amplitude 

Inserting tlie results from (1A.5I) . ( 1A.6I) and fl4.9l) -( l4T2l) into (14. ip . we can write for the 
renormalized one-loop matrix element for a process with n external fermions 

^= MoiQ^) [Ff" + AFr + AFi^] . (4.13) 

Here the corrections are split into a symmetric-electroweak (sew) part, 

•inn \ 2 

='-^EE E /f///(e,Mw;p.,P,)-3C(mOE^o (4-14) 

which is obtained by setting the masses of all gauge bosons, A, Z and W^, equal to Mw 
everywhere, an electromagnetic (em) part 

inn 

= --EE^^^/^ne,0;p.,p,), (4.15) 

^ i=l i=i 

resulting from the mass gap between the W boson and the massless photon, and an 
Mz-dependent part 

= -^EE/f//AJ(e,Mz;p„p,), (4.16) 

^ i=i j=i 

describing the effect that results from the difference between Mw and Mz- The functions 

Al{e,m;pi,pj) = I{e,m;pi,pj) - I{e, Myj;pi,pj), (4.17) 

which are associated with the exchange of gauge bosons of mass m, describe the effect 
resulting from the difference between m and Mw- 

We have expressed the large logarithms in the Yukawa contribution 

C{m,) ""M^ U + \Lle + -Lie' + ^(e^), (4.18) 

which originates only from the counterterm fl4.9p . through their natural scale mt in = 
\n{Q'' /m^), while all the gauge-interaction contributions are written in terms of L = 
ln(Q^/M^) as usual. In order to build squared one-loop expressions to order e°, which 
enter the two- loop amplitude, we need to expand all one- loop terms up to order e^. For 
the functions / and A/ we obtain 

/(6,Mw;p.,p,) -L' - hh - + (3 - 2/,,) (l + hh + h\'^ + 0{e% 

A/(e, Mz;p„p,) ""M^ h {^L + 2Lh + LV) + ^(e^), 

AJ(e, 0; p„ p, ) ^i^^ (5,0 + 5,,o) f-e"' + + + h'eA + (5., + 5,,) (he-' + 



2 6 



2 3 
3 

2^ji - o {^ifl + ^j,o) - 5i,t(l + k) 



e-i + L + iL2e + ^LV)+0(e^). (4.19) 
2 D 
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Only the function AI{e,0;pi,pj), which incorporates the interactions with massless pho- 
tons, depends on the fermion masses rrii through the symbols 6ifi and 5j_t defined in fl2.15l) . 
The functions I{e, My^;pi,pj) and AI{e, Mz;Pi,Pj), which correspond to the exchange of 
massive gauge bosons, are independent of the fermion masses and agree with the results 
for massless fermions in Ref. [32]. 

The functions / and A/ are symmetric with respect to an exchange of the external 
legs and apart from the angular- dependent Zjj-terms, the dependence on pi and pj 
can be separated: 



I{e,m;pi,pj) 



1 r 



I{e,m;pi,pi) + I{e,m;pj,pj] 



(4.20) 



where la = Ijj = is understood. Under the sums over i,j in fl4.13p - fl4.16p we can replace 
I{e,m;pj,pj) — > I{e,m]pi,pi), use the charge-conservation identity fl2.33p and write 



i=l i=i 



i=l 



(4.21) 



and similarly for A/. This relation turns out to be useful later. 



5 Two-loop results 

The renormalized two-loop matrix element gets contributions from the factorizable 
two- loop diagrams (13.41) and (13.151) . from wave- function renormalization, and from pa- 
rameter renormalization, 

M^ = Ml + M^""^ + M^^ + Ml^. (5.1) 

The results for the two-loop diagrams (13. 4p . which involve only gauge interactions, are 
summarized in App. IA.2[ Here we first show that the Yukawa contribution of the fac- 
torizable two-loop diagrams in (13.150 vanishes in NLL accuracy. Then we list the two- 
loop renormalization contributions and finally combine the factorizable contributions from 
two-loop diagrams with gauge interactions and the renormalization contributions into the 
complete renormalized two-loop amplitude. 

As discussed in Sect. [3l the NLL corrections factorize, i.e. they can be expressed 
through correction factors that multiply the Born amplitude. Moreover, as we show, the 
two-loop correction factors can be expressed entirely in terms of one-loop quantities. 
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5.1 Two-loop Yukawa contributions 



As explained in Sect. 13.2.31 the only non-suppressed NLL contributions involving 
Yukawa interactions, apart from counterterms, arise from the first three diagrams in 

(EISD, 




The evaluation of these three diagrams is presented in Sect. IA.3I In NLL approximation 
we find that, up to a minus sign, the integral functions associated with the three individual 
diagrams equal each other, and combining all diagrams we get 

-inn _ I 

i=l j=i Vi=A,Z,W^ *2=i?,X,</>± 

x\Gfur^-ir^Gf^+ y: Gf''i^u}i = o, (5.3) 

where the result for the function Dy is given in flA.531) . Owing to global gauge invariance 
the combination of gauge and Yukawa couplings in the curly brackets of (15.31) vanishes 
[cf. (I2.32p ]. so the only two-loop NLL contributions from Yukawa interactions originate 
from the wave-function counterterms discussed in Sect. 15. 2[ This observation confirms 
the prediction of Refs. [24,25], where in soft-collinear effective theory at scales below Q 
scalar particles contribute only via wave-function renormalization. 

5.2 Two-loop renormalization 

At two loops, the mass renormalization leads to non- suppressed logarithmic terms 
only through the insertion of the one-loop mass counterterms in the one-loop logarithmic 
corrections. However, these contributions are of NNLL order and can thus be neglected 
in NLL approximation. In this approximation also the purely two-loop counterterms that 
are associated with the renormalization of the external-fermion wave functions and the 
couplings, i.e. SZi'^\ Sgf \ Se^"^^ and SXt^\ do not contribute. 

The only NLL two-loop counterterm contributions are those that result from the com- 
bination of the one- loop amplitude with the one- loop counterterms 6Zl^\ Sg\^\ 5e^^^ and 
5\f\ The wave-function counterterms yield 

n 1 

MT = Mli:=^5Z\'\ (5.4) 
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The unrenormalized one-loop amplitude A4j and the wave-function renormalization con- 
stants are given in flA.51) and (14.91) . respectively. In NLL approximation only the LL 
part of Ai^ contributes to (15.41) . At this level of accuracy, the unrenormalized amplitude 
Ail and the renormalized amplitude Ai^ (I4.13P are equal, and we can use (14.211) in order 
to write AiJ and SZ^^^ such that all gauge-group generators If appear only in terms of 
the Casimir operators J2v=a,z,w± ^'^d iflf, which commute with each other. This 
enables us to combine the counterterms Ai"^^ (15.41) with the unrenormalized result At2 
( 1A.46[) into the form presented in Sect. 15.31 

The remaining NLL two-loop counterterms result from the insertion of the one-loop 
coupling- const ant counterterms (14. 4p in the one-loop amplitude (14.131) and read 



^2^PR NLL 



1 

Ye 



Ml, 



i=l 



I{e,Mw;Pi,Pi 



+ 6«e2Q^A/(e,0; 



Pi, Pi 



(5.5) 



Again only the LL parts of the one-loop amplitude contribute to A^2^, so ( 14.2ip has been 
used to arrive at the form (15. 5p . and the Yukawa terms in as well as the corresponding 
counterterms 5X[^^ are irrelevant. 

5.3 Renormalized two- loop amplitude 

Using the results of Apps. |A] and El we find that the renormalized two-loop matrix 
element can be written as 

M2 ^= AioiQ^) I i [Ff + Ff "'AFi""^ + "AFf + ^ [AF^""^]^ + AF^^AFf^" 



+ Gf" + AG 



em I 
2 h 



(5.6) 



in terms of the one-loop correction factors defined in (I4.14I) - (I4.16I) and the additional 
two-loop terms 



J{e,M^^,fj.'^]Pi,Pi), 



1 " 

-y 

1 " f r 

2 E Q' y'-e^ [^Ji^^ 0' f^h p^) - ^Ji^, 0, M^; p„ p, 

+ bi^l^AJ{e,0,M^;p„pM. 



(5.7) 



The two-loop functions J and A J are defined in (1A.48P through the one-loop functions 
/ and A/, so the entire two-loop amplitude is expressed in terms of one-loop quantities. 
The relevant J-functions read explicitly 
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2.^ ^ ^ NLL r / 3 3 , or, -2 , r2,-l 



- (5i,t + 2Lh-^ + 2L^) + 0(e), 

AJ(e, 0,/i^;p,,p,) - AJ(e, 0, M^;p,,pi) "i^"^ j^^.o -2e-^ + e'^l^^ - 2L) 

+ - ^^L] + '^^.t (2Le-i + - \ + 0(e) 



(5. 



where 



Mr 



In 



(5.9) 



Note that the terms ™ and AGl"" in (Q only involve NLLs. 

In order to be able to express (15.61) in terms of the one-loop operators (I4.14I) - (I4.16I) it 
is crucial that terms up to order are included in the latter. 

The coefficients and &qed describe the running of the electromagnetic coupling 
above and below the electroweak scale, respectively. The former receives contributions 
from all charged fermions and bosons, whereas the latter receives contributions only from 
light fermions, i.e. all charged leptons and quarks apart from the top quark. 

The couplings that enter the one- and two-loop correction factorqj are renormalized at 
the general scale /xr > Mw. The renormalization of the coupling constants gi, g2, e, and 
At in the lowest-order matrix element M-q^Q"^) in (14.131) and (15.61) is discussed in Sect. 14. 1[ 
The /iR-dependence of A1q((5^) is imphcitly defined by (14. 8p . and the dependence of the 
one- and two-loop correction factors on /ir is described by the terms (15. 7p . The contribu- 
tions (15. 7p originate from combinations of UV and mass singularities. We observe that the 
term proportional to 6^^-* vanishes for /ir = Mw. Instead, the terms proportional to b^i\ 
62^'' , and cannot be eliminated through an appropriate choice of the renormalization 
scale. This reflects the fact that such two-loop terms do not originate exclusively from 
the running of the couplings in the one-loop amplitude. 

Combining the Born amplitude with the one- and two-loop NLL corrections we can 
write 



where we observe a factorization of the symmetric-electroweak contributions, 

^sew 



NLL ^ _^ «^^sew 
47r ^ 



An) 



2 ri 



- {Fry + G\ 



sew 
2 



the terms resulting from the difference between Mw and Mz, 

pz NLL ^ ^ ^Airz, 
47r 



(5.10) 



(5.11) 



(5.12) 



^ These are the coupHng a in the perturbative expansion (|2.9p and the couphngs gi, g2, e and At 
that appear exphcitly in (|4.14p . (|5.7p and enter imphcitly in (|4.14p - (|4.16p through the dependence of the 
generators (|2.16p on the couphngs and the mixing parameters Cw and Sw 
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and the electromagnetic terms resulting from the mass gap between the photon and the 
W boson, 



pem NLL . Ot 



1 + ^AFr + 



^ '1 {AF^'^f + AG 



'em 



(5.13) 



We also observe that the symmetric-electroweak and electromagnetic terms are consistent 
with the exponentiated expressions 



exp 
exp 



a 



An 



, a \2 

^1 +\^) ^2 



— AFr +(—) ACT 

An ^ \AnJ ^ 



(5.14) 



In particular, these two contributions exponentiate separately. This double-exponentia- 
ting structure is indicated by the ordering of the non-commuting one-loop operators Ff''™ 
and AF^*™ in the interference term /^sew^jT^em result (15. 6p . The commutator of 

these two operators yields a non-vanishing NLL two-loop contribution. 

Note that the 0{a'^) LL contributions in fl5.14p are entirely given by the exponentiation 
of the one- loop term. On the other hand, at the NLL level the presence of the (9(a^) 
terms Gf" and AG|™ in our result (I5.14p seems to spoil exponentiation. However, this is 
an artifact of the fixed-order expansion of the argument of the exponential and must not 
be interpreted as a breaking of exponentiation. This can easily be seen in the framework 
of evolution equations, were the 0{a'^) terms in (I5.14p naturally emerge from the running 
of the coupling associated with the one-loop contribution. To illustrate this feature we 
restrict ourselves to a gauge theory with a simple gauge group and consider a form factor. 
In this case the structure of our result F*^*^™ in (I5.14p is easily obtained from the manifestly 
exponentiated expression (15) in Ref. [30], 



= Fo exp M — 

[J Ml X 



dx 



Ml X- 



— 7(«(x')) + C(«(a;))+e(«(M^)) 



(5.15) 



The NLL approximation requires the one-loop values of the various anomalous dimensions 
as well as the one-loop running of a in 7(a), 



7(a(xO) ^i^^ ^7^^) ^ [1 _ ^^(1) 1, ( 4 
An An An 



7«, 



C(a(x)) ^Ci^\ ^aiM^)) W 

47r 47r 



(5.16) 



Inserting these expressions in fl5.15p one obtains 
^ NLL f «(/^r) [7^^^ r 



^^i^^Foexp/^^ 
[ An 

_ ( <^ ' ' ' 
I An 




(5.17) 
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and one can easily verify that the two-loop term appearing in the argument of the ex- 
ponential, i.e. the term proportional to the /3-function coefficient b^^\ corresponds to the 
term G^^ in our result. 

The one- and two-loop corrections fl4.13p - fl4.16p and fl5.6p - fl5.7p contain various com- 
binations of weak-isospin matrices 1^ , which are in general non-commuting and non- 
diagonal. These matrices have to be applied to the Born amplitude A4q{Q'^) according 
to the definition fl2.3ip . In order to express the results in a form which is more eas- 
ily apphcable to a specific process, it is useful to split the integrals I{e, Mv;Pi,Pj) and 
AI[e,Mv',Pi,Pj) in (14.140 - 04.160 into an angular- dependent part involving logarithms lij 
and an angular-independent part. This permits to eliminate the sum over j for the 
angular-independent parts of fl4.14p - fl4.16p using (14. 21 p . One can easily see that the 
angular-independent part of fl4.14p gives rise to the Casimir operator fl2.18p . After these 
simplifications, all operators that are associated with the angular-independent parts can 
be replaced by the corresponding eigenvalues, and the one- and two-loop results can be 
written as 



where the electromagnetic terms read 

yem NLL ^ ^ !^/\ -)- ( '^'^ | 



em 



(5.18) 



(5.19) 
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(5.20) 



For the term resulting from the difference between Mw and Mz we get 



/ 



Z NLL 



(5.21) 



with 
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Z NLL 
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(5.22) 
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and the symmetric-electroweak contributions yield 



NLL 



a, 



Air 



AttJ 



(5.23) 



with 
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e2 U 
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sew NLL 
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4e2 
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{i)9i 



e2 V2 



(^2 9 



6" 
0(e). 



i=l 



(5.24) 



In the above equations l^^ = ln(/i|^/M^), and Cj, tf, gj represent the eigenvalues of 
the operators Cj, T^, Yi, and Qj, respectively. 

The only matrix-valued expression is the angular-dependent part of the symmetric- 
electroweak contribution /]^^™ in (15.241) . 



i=l j=i V=A,Z,W^ 



(5.25) 



The two-loop corrections involve terms proportional to )Cf^ and [/Cf^J^. However, the 
latter are of NNLL order and thus negligible in NLL approximation. The combination of 
the matrix (15.251) with the Born amplitude, 



-Mo(g^)/cf = EE^^. E < 

i=l j=i V=A,Z,W± 



ipi...(p'-...ip'j...(pri jY jY 



= 1 J = l 



glim, 91.3.3 

e2 4 +e2 ' ^' 



v=w± 



(5.26) 



requires the evaluation of matrix elements involving SU(2)-transformed external fermions 
^'i^^'j, i.e. isospin partners of the fermions fi,(pj. Explicit results for four-particle pro- 
cesses are presented in App. 



6 Discussion and Conclusion 



We have studied the asymptotic high-energy behaviour of virtual electroweak cor- 
rections to arbitrary fermionic processes in the Standard Model. The present analysis 
extends results previously obtained for massless fermion scattering [32] to processes that 
involve also bottom and top quarks. By explicit evaluation of all relevant Feynman di- 
agrams, we have derived a general formula that describes one- and two-loop logarithmic 
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contributions of the form ln(Q^/M^). Such logarithmic terms — which dominate the elec- 
troweak corrections at TeV colhders — originate from ultraviolet and mass (soft/collinear) 
singularities in the asymptotic regime where all kinematical invariants are at an energy 
scale ^ M^. All masses of the heavy particles have been assumed to be of the same 
order Mw ~ Mz ~ Mh ~ rrit but not equal, and all light fermions — including bottom 
quarks — have been treated as massless particles. We have included all leading (LLs) and 
next-to-leading (NLLs) logarithms. 

The calculation has been performed in the complete spontaneously broken electroweak 
Standard Model using the 't Hooft-Feynman gauge. The fermionic wave functions are 
renormalized on shell, and coupling-constant renormalization is performed in the MS 
scheme, but can be generalized easily. Employing the method developed in Ref. [32], we 
have reduced all NLL contributions to factorizable diagrams. In this way the process- 
dependent part of the calculation is isolated in a generic tree-level amplitude, which 
is multiplied by process-independent factors consisting of loop integrals and gauge and 
Yukawa couplings. Technically this is achieved by means of coUinear Ward identities and 
a soft-collinear approximation. All relevant contributions originating from ultraviolet 
singularities are consistently taken into account in every step of the calculation. In par- 
ticular, as discussed in Sect. 13.1.41 the logarithms of ultraviolet origin are isolated in a few 
bare diagrams [see (13.71) and f l3.15p ] and counterterms (see Sects. 14.11 and 15.21) by means 
of minimal subtractions of the ultraviolet singularities and an appropriate choice of the 
subtraction scale. While the ultraviolet logarithms associated with the renormalization 
of the coupling parameters in the factorized tree-level amplitude are easily absorbed into 
running couplings, the process-independent correction factors receive additional contri- 
butions of ultraviolet origin which result from the bare two-loop diagrams of type (13.71) . 
the wave-function renormalization [see (14.121) and (15. 4p ]. and the renormalization of the 
couplings in the soft-collinear corrections [see (15.51) ]. 

Since the NLL electroweak corrections originate only from the electroweak interactions 
of the external legs, our general results for n-fermion processes are also applicable to 
processes involving external fermions and gluons and, more generally, to hard reactions 
that involve n fermions plus an arbitrary number of SU(2)xU(l) singlets as external 
particles. Additional legs associated with external electroweak singlets can only enter the 
hard part (F) of the factorizable diagrams (13. 3p . (13.41) and (13.151) . This modifies only 
the process-dependent hard amplitude M.o, which is always factorized in our derivations, 
while the NLL correction factors receive contributions only from fermionic external legs 
and do not depend on additional external singlets. 

All two-loop integrals have been solved by two independent methods in NLL approx- 
imation. One makes use of sector decomposition, the other uses the strategy of regions. 
Explicit results have been given for all contributing factorizable Feynman diagrams. 

The presence of soft/collinear singularities originating from virtual photons and their 
interplay, at two loops, with logarithmic corrections resulting from massive particles is one 
of the most delicate aspects of the problem. In order to isolate the finite ln{Q'^/M^) terms 
in a meaningful way, one has to separate the photonic divergences in a gauge-invariant 
contribution that can be cancelled against real-photon corrections. To this end, we have 
split the corrections into a finite symmetric-electroweak part, which is constructed by set- 
ting the masses of all gauge bosons equal to Mw, and remaining subtracted parts, which 
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describe the effects resulting from the 7-W and Z-W mass differences. By combining all 
one- and two-loop diagrams we found that these three contributions factorize as described 
in fl5.10l) - fl5.13l) : the term associated with the 7-W mass splitting (F^™) depends only on 
the masses and charges of the external fermions and behaves as a pure QED correction 
subtracted at photon mass Ma = Mw- The term resulting from the Z-W splitting (-F^) 
is proportional to ln(M|/M^) and depends only on the external-leg Z-boson couplings. 
Finally, the contribution constructed by setting Ma = Mi = Mw in all loop diagrams 
^^sew-j ^yj^j^g independent of symmetry-breaking effects such as mixing or cou- 

plings proportional to the vacuum expectation value. This contribution, which contains 
only finite \n{Q'^/M^) terms, behaves as in a symmetric SU(2)xU(l) theory where mass 
singularities are regularized by a common mass parameter. Moreover we find that the 
electromagnetic and the symmetric-electroweak parts exponentiate as described in (15.141) : 
the corresponding two-loop contributions can be written as the second-order terms of 
exponentials of the one-loop contributions plus additional terms that are proportional to 
the one-loop /5-function coefficients. 

These results agree with the resummations that have been proposed in the literature 
and confirm — for fermion scattering processes — the assumption that the asymptotic high- 
energy behaviour of electroweak interactions at two loops can be described by a symmetric 
and unmixed SU(2) x U(l) theory matched with QED at the electroweak scale. Indeed, 
apart from the terms involving ln(M|/M^), in the final result we observe a cancellation 
of all effects associated with symmetry breaking. We have explicitly checked that, upon 
separation of the QED singularities, our results are consistent with the predictions of 
Ref. [19] and Refs. [24,25]. 

As an application of our results for general n-fermion processes, we present in App. [D] 
explicit expressions for the case of four-particle processes involving four fermions or two 
fermions and two gluons. In general, our process-independent results can be applied to 
any reaction with external fermions and gluons as long as all kinematical invariants are 
large. We plan to extend these results to processes involving external gauge bosons and 
Higgs bosons. 

A Loop integrals of factorizable contributions 

In this appendix, we present explicit results for the loop integrals of the one- and 
two-loop factorizable contributions defined in Sect. I3.1.2l and the Yukawa contributions in 
Sect. 13.2.31 These are evaluated within the 't Hooft-Feynman gauge, where the masses 
of the Faddeev-Popov ghosts and would-be Goldstone bosons 0^ read 

M^A = Ma = 0, = M^z = Mz, and M^± = M^w = M^. Using the soft-collinear 
approximation and the projectors introduced in Sects. 13.1.31 and I3.2.2[ we express the 
factorizable contributions resulting from individual diagrams as products of the n-fermion 
Born amplitude with matrix-valued coupling factors and loop integrals. 

The loop integrals associated with the various diagrams are denoted with symbols 
of the type Dh{mi, . . . ,mn;pi,Pj, . . .). The definition of these integrals is provided in 
App. [HI They depend on various internal masses mi,m2, . . . and, through the external 
momenta pi,pj,..., on the kinematical invariants r^j and the masses = pf. The 
symbols are always used to denote generic mass parameters, which can assume the 
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values TJik = My/, Mz, rrit, Mh or mk = 0. Instead we use the symbols to denote non- 
zero masses, i.e. M^. = M-y^/, Mz,mt, M^. The integrals are often singular when certain 
mass parameters tend to zero, and the cases where such parameters are zero or non-zero 
need to be treated separately. We also define subtracted functions 

ADh{mi, . . .,mn;Pi, • • •) = Dh{mi, . . .,mn;Pi, ■ ■ ■) - Dh{Mw, ■ ■ .,Myi/;pi, . . .), (A.l) 

where the integral with all internal mass parameters equal to Mw is subtracted. 

The integrals have been computed in NLL accuracy, and the result is expanded in e 
up to C(e^) at one loop and C(e°) at two loops. The UV poles have been eliminated by 
means of a minimal subtraction as explained in Sect. 13. 1.41 such that the presented results 
are UV finite. The integrals have been evaluated separately for all physical combinations 
of gauge-boson and fermion masses on internal and external lines. All loop integrals have 
been solved and cross-checked using two independent methods: an automatized algorithm 
based on the sector-decomposition technique [33] and the method of expansion by regions 
combined with Mellin-Barnes representations (see Ref. [31] and references therein). 

The one-loop diagrams are treated in Sect. IA.lt the two-loop diagrams involving gauge 
interactions in Sect. IA.21 and the diagrams involving Yukawa interactions in Sect. IA.3I 

A.l One-loop diagrams 

The one-loop factorizable contributions (13. 3p originate only from one type of diagramj^ 




MY= IXfC = -^0 E irirDoiMy,;p.,,p,). (A.2) 

Vi=A,Z,W± 

~3 

In NLL accuracy, the representations of the generators 1^^ and ij^ correspond to the 
chiralities given by the spinors of the external particles i and j, respectively. The loop 
integral Dq is defined in f lB.6|) and to NLL accuracy yields 



NLL j2 2 3 1 r4^2 , (a oI \ f T , ^T^^ , ^r3 2 



Do{Mi;p„p,) "^"^ -L' - -L'e - -L\' + (4 - 2/,,) L + -L't + -L't 

3 4 V 2 



+ h [21 + 2Lh + 
DoiO■,p^,p,) 2k,e-' - (5,,o + 5,,o) (e'' + 2e-') + I 5,, 



2 6 



+ iL'e' - he-' + {2-k)(L + -Lh + -Lh' 
24 ^ ' \ 2 6 



+ (z ^ j) 



(A.3) 

where the UV singularities 

D^''{rm;p,,p,)''M^4e-' (A.4) 



^The Z-loop diagrams depicted in this appendix are understood without factors (ae/47r)' 
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have been subtracted. The shorthands L, / j, kj, 5j_t, 5i,o are defined in Sect. 12.11 

Summing over all external legs, we find for the factorizable one-loop contributions 
(1^ 



with 



-inn 

^ i=i j=i v=A,zyv± 

1 n n 



FFjSew 
1 

2 

^ i=i j=i 

^^1^'" = -2EE^f^/A^o(Mz;p„p,- 



5:5:/f//ADo(0;p.,p,), 



and ADQ{m;pi,pj) defined in (lA.ip . 

A. 2 Two-loop diagrams involving gauge interactions 



(A.5) 



(A.6) 



The two-loop NLL factorizable terms (13.41) involve fourteen different types of diagrams. 
The diagrams 1-3, 12 and 14 in this section give rise to LLs and NLLs, whereas all other 
diagrams yield only NLLs. 

Diagram 1 




NLL 



M, J2 lf'lf'lJ'lJ'D,{Mv,,Mv,;PuPj), (A.7) 

Vi,V2=A,Z,W^ 



where the loop integral Di is defined in (IB.6P and yields 

Di{Mi,m2;pi,Pj) 
Di{0,M2;pi,Pj) (4o + M 



6 3 ^ u ^ ' 



4 

3^ 



(4 - 2/,,) ( + L^e-i + 



-Lh-^ - + (4 - 3/i, + 2h + /i) i^'e-i 



'20 17, 16, 7, ^ . 
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3^-1 



9 



-L'e 



36 



4 _ 4 2 \ _3 
3 3'^^ + 3'^'' 



(4 - - /.) [he-' + + 



+ (i ^ j) 



2L\-' + ^Lh-' + + 4/,,Le-2 + (8 + 6/.,) L'e~' 
3 2 ' -^^ 



Here the UV singularities 



NLL at2-1 



have been subtracted. 



(A.9) 



Diagram 2 



Ml''' 




NLL 



M 



E lf^lfujHfD^{Mv,,My,-p,,p^), (A.IO) 



where the loop integral D2 is defined in (1B.6I) . This integral is free of UV singularities 
and yields 



02(0, M2; Pi, Pj) (5i,o 



^lV^ - + (4 - 2y (lV^ + l^ 



+ I2 { 2L'e-' + 



h^-[4-^l^^ + 2l2 + h,]L' 



D2{MuO;pi,pj) = D2{0,Mi;p,,p 



D2{0,0■,p^,p,)''= 6,,o6,,o e-' + {4-2h,)e-'' + <! M.,o 



_ 1^6-3 + lL2e-2 
6 3 3 



+ -L'e-' + — + (2 - 2/i,- + k) -e 



-Lt 



-2 



4 2. 



3 3 



+ ( :^ + ■::k, - ik ) i^^e"^ + ( ^ + -k, - ^k) L 



16 2. 



10, 



9 9 
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Diagram 3 



3 3*^ 



+ (/i + /,)(2LVi + ^L3 



3 3^ 



(A.ll) 




NLL •fl'2 . . 



^ ey^v^v^lfufufWs{Mv„Mv„Mv,;p^,Pj), (A.12) 
Vi,F2,v'3=A,z,iy± 



where the e-tensor is defined in Ref. [32] [see also fl2.17p ]. The loop integral is defined 
in (1B.6I) and yields 



D3{Mi,m2,M3;pi,pj) ^ 
Ds{0,M2,Ms;p^,Pj) ""=6, 



^^L' + (2 - + h) Lh-' + (1 - 
3 V3 3 '^^ 3 7 



D3{M,,M2,0;p„p,) "i^^ -35,,o^e-2 + 6,,, (^^'e"^ + ^ 



12 



+ 5 



where the UV singularities 



3Le-2 _ (2 + - I,) Lh-' + - + ^/i) 

+ + + -h- ll,) L'] , (A.13) 



D^^(mi,m2,0;pi,pj) '^=^ -3 (5i,o + ^j.o) e"^ 

have been subtracted. 
Diagram 4 



(A. 14) 




v.v^.^A.,H^- (A. 15) 
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where the loop integral D4 is defined in ( ]B.6|) and yields 



-S,J-L\-' + -L' 
^' V2 6 



( 2Le-2 + sLh'' + yL^ ) . (A. 16) 



Here the UV singularities 



l^r(0, =^ (5,,o + 5„o) - (5,, + 5,,t) (^e-^ + Uh^' + ^L^) (A.17) 

have been subtracted. 

Diagram 4 and the following diagram 5 are the only cases where fermion masses in the 
numerator of the fermion line i contribute to the result in NLL accuracy. In principle, these 
fermion-mass terms contribute with generators 1^ and 1^ , belonging to representations 
with different chiralities [see f l2.29l) ]. But we found that this happens only in the case when 
both gauge bosons Vi and V2 are photons, such that if = if. Thus all contributions to 
(IA.15I) and flA.181) can be expressed in terms of the operators Ij , Ij, which belong to the 
representations associated with the chiralities of the external fermions. 

Diagram 5 




Mr= i®cr^ = -A^o E iriririrD,{My,,My,-p,,p,), 

v.y.^A,z,w^ (A.18) 

where the loop integral is defined in ( IB. 61) and to NLL accuracy is given by D4, up to 
a minus sign: 

Dr,{mi,m2;Pi,Pj) -Di{mi,m2;pi,Pj). (A.19) 

Note that this relation only holds if the fermion-mass terms in the numerator along the 
line i are correctly taken into account. 
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Diagrams 6 




+ 




^1,^2^3,14=^,^,1^* (A 20) 



where the loop integral is defined in f lB.6P and yields 

DQ{Mi,m2,m3,M4;pi,pj) -^L^, 

y 



D6(0,M2,M3,M4;Pi,p,) 



NLL 20 3 M| + M| 



9 



2M| 



D6(Mi,M2,M3,0;pi,p,) = 1^6(0, M2,M3, Mi; 



1^6(0, M2,M3,0;p„p,) "i^^ (5^,0 + M (^yLe'^ + ^LV^ 

Here the UV singularities 

Dr(Mi,m2,m3,M4;p.,j>,) ^i^^ ^^V^ + 

NLL 10 2 -1 , 20 3 ml + ml 



10 



(A.21) 



/^^^(O, ms, mg, M4; Pi, pj) 



3 9 



2MI 



{Ml, 1712, 771^,0; Pi, pj) = ^6 {0,rn2,m3,Mi;pi,pj) 
D^'^{0,m2,ms,0;pi,pj) y (5i,o + 5i,o) 



(A.22) 



have been subtracted. We observe that the loop integrals associated with A-Z mixing- 
energy subdiagrams give rise to the contributions 



ADe{0,M^,M^,Mz■,p^,PJ: 



NLL 



(A.23) 



which depend linearly on the ratio M^/M^. Similar terms appear also in diagrams 7, 8, 
9, and 10. These terms cancel when adding all contributions owing to relations between 
these integrals (see App. [C]), which hold also in the presence of massive external fermions. 
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Diagram 7 




E 



i. j 



^ViVzV^ViViV 



Vi,V2,V3=A,Z,VK± 



V=A,Z,W^ 



x{D- l)Dr{Mv„Mv,,Mv,;Pi,Pj), 
where D = A — 2e. The loop integral Df is defined in (IB.6P and yields 



Dr{Mum2,M3■,p^,PJ) 0, 



Dr{MuM2,0;p„pj) = Dr{0, M2, M^;pi,pj), 
Dj{0,M2,0;p,,Pj) ''i^^O, 
where the UV singularities 



D}/''{M,,M2,0;p^,Pj) = Dr(0,M2,Mi;p,,p,) 
have been subtracted. 



(A.24) 



(A.25) 



(A.26) 



Diagram 8 




X 



{/^«, J^^}^ ^ D8(Mv^,, Mvi; 



(A.27) 



where the curly brackets denote anticommutators and v is the vacuum expectation value. 
The loop integral is defined in (IB.6P and yields 



M:^Ds{M,,M2,m^,M^;p,,p,) ^i^^ 0, 
Ds{0,M2,Ms,M,;p,,p,)''M^ ^ 



Ds{M^,M2,M^,0;pi,p,) = Ds{0, M2, M3, M^;pi,pj), 



M^Ds{0,M2,Ms,0;p„p,) 0. 



(A.28) 
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Here the UV singularities 



15^(0, M2, Ms, M4;p„p,) ""^"^ 



have been subtracted. 



(A.29) 



Diagram 9 



Ml''' 



^ V . > NLL 1 . . 

:JfI *3< = Mq 

'^i 2 Vl,V4=A,Z,W± 



E 



Vi 



« 3 



X DQ{Mv,,M^2,M^^,Mv^;Pi,Pj), 
where the loop integral Dg is defined in (IB .6 1) and yields 



(A.30) 



Dg{M^,M^,M^,M,-p,,p,) ^i^^ 



M2 + M2 



9 ' 2M| 



- (5,,t + 5„ t) (2L^e"i + 3L 
A(Mi,M2,M3,0;pi,p,) = D9(0,M2,M3,Mi;pi,p,), 



NLL , r N / 1 r -2 ,1 r2,-l 



(A.31) 



Here the UV singularities 



DUV(Mi,M2,M3,M4;p„p,) 



DUV(0,M2,M3,M4;p,,p,) 



3 9 ' 



9 



2M| 



'2,-1 



Dr(Mi,M2,M3,0;p,,p,) = D9^^(0,M2,M3,Mi;p,,p,) 
DUV(o, M2, M3, 0;p,,p,) "i^^ ^ (5,,o + ho) e"' 



(A.32) 



have been subtracted. 
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Diagram 10 




NLL 



Vi,V3=A,Z,W± 



4>2*2 



where 



Diagram 11 



D 



10 



Dr. 



(A.33) 
(A.34) 



For the diagrams involving fermionic self-energy subdiagrams we consider the contri- 
butions of a generic fermionic doublet \1/ with components = u, d. The sum over the 
three generations of leptons and quarks is denoted by Y^^n and colour factors are implic- 
itly understood. Assuming that all down-type fermions are massless, = 0, and that 
the masses of up-type fermions are = or mt, we have 








E 



i j 



Vi,Vi=A,Z,W^ 



X 



E^ E ^ ll^-^-ll"-^- Dii^Q{Mv^,mi^,mi.^,Mv^]Pi,Pj) 
llk^nlX^+lXjXn)mlD^^^,^{My,,mu,^ (A.35) 



where Du m = — 4-D8 represents the contribution associated with the m„-terms in the 
numerator of the up-type fermion propagators of the loop insertion, whereas the integral 
Dii^Q, which is defined in (1B.6I) . accounts for the remaining contributions. This latter 
integral yields 

Diifi{Mi,m2,m3,M4;pi,pj) ^L^, 
Dn,o(0, m,, ms, M,;p,,p,) ^i^^ + ^^i±^ + 5.,o) {^Le-' + 21^-' - -^L') 

- (5i,t + 5,,t) {^Lh-^ + QL^ 
Du^o{Mi,m2,m3,0;pi,pj) = Du^o{0,m2,m3, Mi;pi,pj), 
Dn,o{0,M2,M,,0;p„Pj) ""i^^ (5.,o + M (^Le'^ + ^lV^ 

Dn,o(0, 0, 0, 0; p„ pj) ""i^^ - ((5,,o + 5,,o) e"' + i^t + <5j,t) (^i^f-^ _ , (A.36) 
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where the UV singularities 



2,-1 



D^;^ Q{Mi,m2,ms,0;pi,pj) = D^^ Q{0,m2,ms, Mi;pi,pj) 
D^^q{0, 1712, 1713,0; Pi, pj) ^{5i,o + 5jfl)e~^ 



(A.37) 



have been subtracted. As a consequence of 



2 I 2 

ADu^o{0, m2, ins, M4; Pi, pj) — ^ AL'ii_^(0, mg, mg, M4; p^, p^), 



AL)ii,o(Mi, m2, m3,0;pi,pj) 



NLL "^2 + 



ADn^rn{Mi, m2, m^, 0]pi,pj), (A.38) 



all terms proportional to the fermion masses in A^g^'*"* cancel. 
Diagram 12 




M, J2 If'I^'IJ'Il'D,2{Mv,,Mv,■,p^,p„Pk) 

Vi,V2= A,Z,W± 



(A.39) 



where the loop integral D12 is defined in (IB.6P and yields 



D,2iM,,M2;p^,p„Pk) - (4 - 2hk + ii + ^/a) i^', 

12 



Dl2(0, M2; Pi, Pj, Pfc) ^= {^ifi + (^j.o 



- (2 - {2Le-^ + L^e-i + 11^) - ^2 (2Le-2 + ^.Lh'^ + 



-^L^'e-i - ^-L" + (2 - 2Z,fc + + /.) L\-^ + {^-- 2kk + ^^2 + ) 



+ 6 



4 



3^-1 



7 



+ (6 - 2/ifc - 2/,j + 3/2 + Ij) L'e 



2,-1 



36 



Di2{Mi,0;pi,pj,pk) Sk,o 



-U'e-' - + (2 - k,) ( 2L\-' + 



3 

—L^ — ( lik H — li H — lk \ L"^ 

6 V 3 3 3 3 7 



2e-^ + 



2 4 2 2 1 4 

3 3 3 9 3 V jfc fc; 



- (4 - - 4) [he-' + ^L^6-^ + 



_ + ^L^^e-^ + -^L^ + (2 - 2/^^ + k) e'^ - (4 - 2/^^) Le" 



+ {kk - k) L'e-^ + ( 3 + 3^*-^ -U]L' 



4 1 



+ 4,t 



1 1 1 

3 12 



1, 1 



4) - (4 - 2kk) 



-L'e 



+ 4,1 



3 



_1 1 5 , _3 

Hk 3 nj ~l~ g J 6 



4 — -/ifc — -kj + -Ij — -Ik ] Le ^ 

,'4 7, 10, 5, 4 . 3 
+ [---lik + Y^ij-Qh-Q^k]L 



(hk - ^hj + 21 J + Ik) L'e 



2,-1 



Le"^ + L^e-' + 21^6'' + 2L*- { I, 



3^-1 



2' 2' 



(4 - 4kk + 2k) Le-' + (4 + 2kk + 2/,, - 5/, - /,) L^e 



2-1 



10, 



19 



+ 4,1 



2L'e-' + 2L^e-^ + -L^ + {2kk + 2/^,- - 2k - Ij - h) Le 



+ (8 - 2kk + 4/i,- - 2>k - 2lj - Ik) L'e-^ + 8 - Akk + —kj - i:k - i^lj L 



14, 



(A.40) 
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Here the UV singularities 



(Ml, 1712] Pi, Pj,Pk] 



D^^ (0,1712; Pi, Pj,Pk] 



NLL 



NLL 



-4 (5i,o + 5j,o) 



_-3 



(A.41) 



have been subtracted. While the above diagram, to NLL accuracy, does not depend on 
Tij and Tjk for p| = 0, a dependency on r^j is introduced for p'j = ml- 

Diagram 13 




NLL . 92 







E 



e 



V1V2V3 jVi tV2 jVi 

-'j -'fc 



Vi, V2,V3= A, Z,W± 



X Dis {Mv, , , ;Pi,Pj,Pk), 



(A.42) 



where the loop integral D13 is defined in (IB.6p . This integral is free of UV singularities 
and yields 



Dn{M^,M2,M,;p„p„pk) ""=0 



Dis{0,M2,M3;pi,pj,pk) {lij - I 



ik 



+ ^kt il2-h)L\ 
D^s{Mi,0,M^;Pi,Pj,Pk) = D^^{0,M^,Mi;p„pk,Pi) 
Di3{M,,M2,0;pi,Pj,Pk) = Di3{0,Mi,M2;pk,Pi,Pj) 



(A.43) 



When one of the gauge bosons is a photon, it couples to two W bosons with equal masses, 
so the terms with the mass- dependent logarithms li, I2, 13 in flA.43p vanish in all physically 
relevant cases. 

Diagram I4 



Ml"'''"' 




NLL 



M, lf'lJ'lluY'D,,{My, , My, ; p„ p„Pk,Pi) 

Vi,V2=A,Z,W± 



{AAA) 

where the loop integral D14 is simply given by the product of one-loop integrals flA.3p . 

Du{Mv^, Mv,;Pi,Pj,Pk,Pi) = Do{Mv^;Pi,Pj) Do{Mv,;Pk,Pi)- (A.45) 
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Sum of two-loop diagrams involving gauge interactions 



The complete contribution of all factorizable diagrams not involving Yukawa contribu- 
tions is obtained by inserting the above results into (13 ■4p . For the case of massless external 
fermions, we have explained in detail in App. E of Ref. [32] how the factorizable two-loop 
diagrams can be summed up to the total two-loop amplitude. To this purpose we have 
used relations between the scalar loop integrals which are listed in App. B of Ref. [32] 
and are valid in NLL approximation. For the general case of diagrams with massive and 
massless fermions, these relations receive only minor modifications which we indicate in 
App. Oof the present paper. Apart from that, the whole procedure remains exactly the 
same. So here we only present the result and refer to Ref. [32] for more details. 

The factorizable two-loop contributions can be written in the form 



M 



F NLL 



1 r 



F,sew 



^F,em 



F,scw 



iF,Z 



+ 2 



af; 



F,em 



nF,em 



F,sew 



iF,em 



(A.46) 



where the one-loop terms are given in ( lA.6p . The additional two-loop terms read 



2/^F,sew 



1 " 

-y 



b?g![f] +b',"gia 



J{e,My^,Q'^]Pi,Pi 



AG 



F,em 



J2 Q\ b^e^ fA^(e, 0, Q''■,p^,p^) - A J(e, 0, M^; p,, p., 



i=l 

'(1) A T^, n i\^2 



+ b\^>,^AJ{e,0,M4,;p„p,)\, 



(A.47) 



with the one-loop /5-function coefficients (14. 5 p and (14.61) . the SU(2) Casimir operator 
(I2.18p . and the two-loop functions 



J{e,m,^'^;Pi,Pj] 



I{2e,m;pi,pj) - ( — 1 I{e,m;pi,pj) 



AJ{e,m,iJ,'^;Pi,Pj) = J{e,m, fi"^; pi,pj) - J{e, , n^;pi,pj), 



(A.48) 



for m = Mw,Mz,0, which are combinations of one-loop functions /. The expressions 
flA.47p rely on the fact that the J-function, to NLL accuracy, involves only the LL parts 
of the /-function. In particular, no angular-dependent /jj-terms are relevant for the J- 
function, so the identity (I4.20p yields 



J{e,m,n'^;Pi,Pj 



NLL 



J{e,m,fi^;Pi,Pi) + J{e,m, fi^;pj,pj) 



(A.49) 



and ( 14.2ip can be generalized to the functions J and A J. 
In order to combine the terms in ( 1A.47P with (15. 5p we use 



-VI -1 



I{e,m;pi,pj) = J{e,m, ij,^;pi,pj) - J{e,m,Q ■,Pi,Pj) (A.50) 



and a corresponding relation between AI and AJ. 
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A. 3 Yukawa diagrams 

Most diagrams witli scalar bosons coupling to external fermions are suppressed, as 
explained in Sect. I3.2.3I The only relevant Yukawa contributions from bare diagrams are 
presented in the following. 

Yukawa diagram 1 




1 



X DY^i{Mv,,M^,;Pi,Pj). (A.51) 
The loop integral Dy i is defined in flB.6p and yields 

DY^i{mi,M2;Pi,Pj) DY{mi;pi,pj) (A.52) 

with 

- {^^Lh-^ + ^^^) , (A.53) 
where the UV singularities 

Dy {Mi;pi,pj) = l^L e + 3^ ' 



(A.54) 



have been subtracted. 
Yukawa diagram 2 




X DY,2{Mv„M^,;Pi,Pj). (A.55) 
The loop integral -Dy,2 is defined in (]B.6|) and yields 

DY,2{mi,M2;Pi,Pj) -DY{mi;pi,pj) (A.56) 
with DY{rni]pi,pj) from (1A.53I) . 
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Yukawa diagram 3 




X DY,z{M^,.M^,,Mv,]Pi,Pj). (A.57) 



The loop integral -Dy,3 is defined in flB.6P and yields 



DY^s{Mi,M2,m3;pi,pj) -DY{rn^;Pi,Pj) (A.58) 
with DY{ms;pi,Pj) from flA.53p . 

B Definition of the loop integrals 

In this appendix, we list the explicit expressions for the Feynman integrals that con- 
tribute to the one- and two-loop diagrams discussed in App. [A] In order to keep our 
expressions as compact as possible we define the momenta 

ki=Pi + li, k2=Pi + l2, k3=Pi + li + l2, 

qi = pj - k, q2 = Pj - k, qs = Pj - h- h, h = -h - h, 
ri=Pk-k, r2=Pk-l2, rs=pk-li + l2, k = h - k- (B.l) 

For propagators with mass m we use the notation 

P{q, m) = q^ -m^ + iO (B.2) 

and for triple gauge-boson couplings we write 

r>^^^'>^'{h, k, h) = g'^'^'ih - hr + g^'^'ih - + g^'^'ih - hT- (b.3) 

The normalization factors occurring in (12. 9p are absorbed into the integration measure 
and for the projection introduced in fl3.12p we use the equality 



n., (^.T) = in,, (F) , (B.5) 

which holds if F does not involve 75 or cur^l, as it is the case in the following equations. 
We have explicitly verified that the NLL contributions from the projector Ilj, in (13.131) 
cancel. 

The integral functions Dh depend on the internal masses mi, 1712, ■ ■ ■ and, through the 
momenta pi,pj, . . ., on the kinematical invariants Vij and on the masses mf = pj of the 
external particles. The definition of these integrals also involves the masses m^, mf , m'-' of 
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the particles </3^, ^p'l, ^p'l' along the fermionic line i after one, two or three interactions with 
gauge bosons or scalar bosons, which possibly change the weak isospin of the external 
particle </?j. We have found, however, that the dependence of the results on the masses 
along the fermionic lines is completely fixed by the external masses m,, so we do not 
indicate the additional masses m[, ... in the arguments of the functions D^. 
With this notation we have 

Aikiqi 



Do{mi;pi,pj)^ dh— , 
J P(/i,mi)P(fci,m;)P(gi,m;.) 

Di{mi,m2;pi,Pj) = / dlidk 
D2{mi,m2;pi,Pj) = j dhdk 



' P{h, mi)P(/2, m2)P{h, m',)P{h, m'l)P{q,, m'j)P{q2, m]) ' 

-16(/i:ig3)(A:3g2) 

' P(/i, mi)P(/2, m2)P{h, m'^P{h, m'l)P{q2, m'^)P{qs, m'j) ' 
D3{mi, 7112, 7713; Pi, Pj) = 

/ .J .7 -2n,,- [(1^3 + m/Dr^ip, + m^T^I q^'^^^^..,^s{h, k, k) 
J ' ^P(h, m{)P{k, m2)P{k, m3)P(h, m'^P(ks, )P(g3, m'j) ' 

D.im.,m2;p.,p,) ^ I dUk ^n, + rnDr^l^s + <)7..(l^. + O^:] 
41 2,Pr,Pj) J 'p{h,m,)P(k,m2)[P(h,m'^fP{h,m^)P{q„m'j)' 

u,[m,,7n2, p^,pj) J ai,ai2 ^^^^^ ^^)p^i^^ m2)P{h, m'nP{h, m'l)P{k2, m[)P{q„ m^) ' 
DG{mi,m2,m3,7ni;pi,pj) = 

-4A:rgiM4 [^^^,.^Ah, k, k)T^'^'^'-^{h, k, k) + 2/2^1/3^*] 



J dhdk 



P{h, mi)P{k, m2)P{k, m3)P{h, mi)P{ki, m[)P{q^, m'^) ' 



D'j{mi, 1712,7713, Pi, pj) = j dlidk 



''P{h, mi)P{k, m2)P(/i, m3)P(A;i, mOP(gi, m'^) ' 
Ds{7ni, 7712, 7713,7714; Pi, Pj) = 

-Akiqi 



J dhdl2- 



■P{h, mi)P{k, m2)P{k, m3)P{h, m^jPiki, 77i[)P{qr, 7n'j) ' 
Dg{mi,m2,m3,m4;pi,Pj) = 

.7 .7 ^k^i'q^\k-kUik-kU 

' ^P{h, 7n,)P{k, m2)P{k, m3)P{h, m4)P{k,, 7n'i)P{qi, m'^) ' 

Dio{mi,m2,m3;pi,pj) = D7{mi,m2,m3;pi,pj), 

Dii,o{mi,m2,m3,m4;pi,pj) = 

4A;f^gr Tr (7^ J 27^4/3) 



J dfid/2 



P{h, 7ni)P{k, m2)P{k, m3)P{h, m4)P{ki, 7n'i)P{q^, m'j) ' 
Dn,m{mi,7n2,7n3,7n4;pi,Pj) = -ADs{mi,7n2,m3,7nA;pi,Pj), 
Di2{mi,7n2;pi,Pj,Pk) = 

-16(/i;igi)(A;3r2) 



j dhdk 



P{h, mi)P{k, m2)Piki, 7n'i)P{k3, 7n'l)P{qr, 7n' )P{r2, m'J ' 
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' ir iT 8A;f^g^V^T^^^3^3(-/i, I2, h) 

' ' P{h, m,)P{l2, m2)P{h, m3)P{h, m'^P{q2, m'^)P{rs, m',) ' 
Du{mi,m2;Pi,Pj,Pk,Pi) = Do{mi;pi,pj)Do{m2]Pk,Pi), 

D (ra m-v v ) - f dl dl -^^^Aih + m^ih + ^'Dih + mdii] 

^''^ ' " ^ ' 'P{h,rm)Pil2,m2)[Pih,m'^fP{h,m1)P{q,,m'^) 

DY,2{mi,m2]Pi,Pj) = 

f .7 .7 -2% [(^1 + m'nih + m'l)m2 + mQ] 

Dy^S^mi, 1712,1713; Pi, Pj) = 



Jdk 



,r 2n,,P3 + O(^i+m0]g3/4 . . 

^P{h,mi)P{l2,m2)P{h,ms)P{h,m'^P{h,m'{)P{qs,m'^)- ^ ' ^ 



The previous definitions are valid for diagrams witli incoming fermions. In tlie case of 
an incoming antifermion (pi, all masses mi,m[, . . . along the fermionic line i have to be 
multiplied by (—1), in the integral definitions flB.6P as well as in the projectors (13.121) . But 



as mentioned in Sects. I2.2l and r3.2.2[ the NLL results are insensitive to this transformation. 

The various coupling matrices and Gf , which are associated with the interactions 
along the fermionic lines k = have been factorized from the loop integrals and 

can be found in App. |Al Here a comment is in order since, in principle, the fermion-mass 
terms in the numerator flip the chirality of the fermions and give rise to coupling matrices 
l]^ and Gf corresponding to opposite chirality states [see fl2.29p ]. However, as discussed 
in App. El for the case of diagrams 4 and 5 [see text after (1A.17I) ]. we have found that the 
fermion-mass terms in the numerator are only relevant in the case of photon interactions, 
where the representations of the generators are independent of the chiralities, i.e. = I^. 
Thus, all contributions can be expressed in terms of the operators which belong to 
the representations associated with the chiralities Kk of the external fermions, and all 
coupling factors can be factorized as in App. \^ 



C Relations between loop integrals in NLL approximation 

In order to combine the two-loop contributions of Sect. IA.2I we use relations between 
the loop integrals. These relations have been obtained from the explicit results listed in 
Sects. lA.ll and IA.2I They are valid after subtraction of the UV singularities and in NLL 
approximation. 

For the case of massless fermionic particles, the relevant relations have been listed in 
App. B of Ref. [32]. We have found that, after only small modifications, they are all still 
valid for the case of massive fermionic particles. One trivial and obvious modification is 
the following change of arguments in the integral functions Dh. 

Dhi- . .■,rij) ^ Dh{. . .■,Pi,Pj), ioi Dl, . . . , Dio, Dufl, Du,m, (C.l) 

and similarly for the subtracted functions AD^ defined in (lA.ip . Many relations do not 
need further modifications, and we refer to App. B of Ref. [32] for them. 
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However, since the presence of fermion masses breaks the invariance of some diagrams 
with respect to an exchange of external or internal lines, certain relations obtained for 
massless fermions have to be modified by an appropriate reordering of arguments in 
the D/i-functions. We list these relations in the following. As in App. \^ the symbols 
rrii are used to denote generic mass parameters, which can assume the values rrii = 
Mw, Mz, rrit, Mh or = 0, and the symbols Mj are used to denote non-zero masses, i.e. 
Mi = Mw,Mz,mt,MH. 

In the second line of (B.2) in Ref. [32] the arguments Pi,Pj have to be exchanged on 
the right-hand side. This relation becomes 

D2{mi,m2]Pi,Pj) = D2{m2,mi;pj,pi). (C.2) 

In the relations (B.3) of Ref. [32] the order of the mass parameters in D2 has to be 
reversed, and the order of the momenta in the 3- leg integral Du is now important: 

Dsi^Mi, 1712, 1713; Pi, pj) ^D2{Mi,m3;pi,pj) - 04(1713, Mi, pi,pj) 

- 6Dg{m3,Mi,Mi,m3;pi,pj), 
AD3{Mw,Mw,mi;pi,pj) -/\Di2{M^,mi;pi,pk,Pj) - AD^^rrii, MY^;pi,pj) 



- 6ADg{mi,M^^,My^,mi;pi,pj), 



AD3{mi, M^ , M^;Pi,Pj) ^= AD3{My^,mi, My^;pi,pj) + ADi{My^,mi;pi,pj) 

+ AD2{mi,Mw;Pi,Pj) + AD4{Mw,mi;pi,pj) 

- -ADi2{Mw,mi;pj,pk,Pi). (C.3) 

As in Ref. [32], the first of these relations has been verified and is needed only if at most 
one of the masses m2 and 7123 is zero. 

The functions J and A J defined in (IA.48P of this paper now also depend on the 
external momenta Pi,Pj, and (B.6) and (B.7) of Ref. [32] become 

3DsiMi,M2,M3,M4;Pi,Pj) -J{t,M^,Q''-pi,pj), 
3AD,{0,M2,M3,0;pi,Pj) - [AJie,0,Q^;pi,pj) - AJie,0, M^;pi,p,] 

3 [ADu,o{0, 0, 0, 0; pi, pj) - ADn,o(0, M2, M3, 0; pi, pj)] -4A J(e, 0, M^; pi,pj). 

(C.4) 

For the 3-leg integrals the order of the external momenta pi,pj,pk becomes relevant 
and (B.8), (B.9) and (B.IO) of Ref. [32] generalize to 

Di2{mi,m2;pi,Pj,Pk) + Di2{m2,7ni;pi,pk,Pj) Do{mi; pi, pj) 00(1712; Pi, Pk) , 

J2 sgn{7r{z, J, k)) Dum,M2;p.,p„Pk) 0, (C.5) 

TT{i,j,k) 

where the sum runs over all permutations 7^(1, j, k) of k, with sign sgn(7r(z, j, k)), and 
Di3{MuM2,M3;pi,pj,Pk) 0, 
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2 ADis {Ml, Ml, ms; Pi, pj , pk) 



2AL)i3(m3, Mi,Mi;pk,Pi,Pj) 
2ADi3{Mi,m3, Mi;pj,pk,Pi) 

ADi2{M^^,m^;Pj,Pi,Pk) - ADi2{M^^,m^;pi,pj,pk). 



NLL 



(C.6) 



When fermion masses are involved, the last relation is only true (and only needed) if the 
two non-zero gauge-boson masses Mi on the left-hand side are equal. 

We have found that the relations from App. B of Ref. [32], together with the modifica- 
tions presented above, are exactly the ones needed to combine the two-loop contributions 
of Sect. lA. 21 into the complete amplitude (see end of App. 

D Application to four-particle processes 

Here we apply our results to four-particle processes involving light fermions, heavy 
fermions and gluons. We first examine four-fermion processes 



where each of the tpi may be a fermion, ipi = or antifermion, (pi = with the 
notations from Sect. [21 provided that the number of fermions and antifermions in the 
initial and final state is equal. We exclude top quarks from the initial state, ipi^2 7^ t, i, 
but allow for bottom quarks there. The final state may contain any combination of 
massless or massive fermions and antifermions, including top and bottom quarks. 

The scattering amplitudes for the processes (ID.ip follow directly from our results for 
the generic n — process fl2.ll) by crossing symmetry, and the Mandelstam invariants are 
given by s = ri2 = T34, t = ri3 = r24, and u = = with rij = {pi + pj)"^. In practice 
we can restrict ourselves to the calculation of s-channel amplitudes f^^ /^^ ~^ fal fat^ 
amplitudes where the external fermion lines are connected between f^^ and /^^ 
initial state and between f^^ and fl^^ in the final state. All other scattering amplitudes 
needed for the four-fermion processes can be obtained from the s-channel amplitudes by 
crossing symmetry. 

In Sect. ID.ll we treat neutral-current four-fermion amplitudes where the particle pairs 
in the initial and final state are antiparticles of each other. Sect. ID.2l is devoted to charged- 
current four-fermion amplitudes where the initial and final state each consist of a pair of 
isospin partners. Finally, in Sect. ID.3l we provide results for the annihilation of two gluons 
into a fermion pair, gg — >■ fa- -^^^ other four-particle processes involving two gluons 
and two (ant i) fermions, i.e. gf^ g/^) fa fa ^ S&y ^tc, are related to gg — > /^/^ by 
crossing symmetry. 

In order to keep all results manifestly invariant with respect to crossing symmetry, we 
keep the hard scale Q^, which enters the logarithms L = \n{Q'^/M^), as a free parameter. 
In practical applications, can be identified with the centre-of-mass energy s or, alter- 
natively, with |t| or \u\. This implies an ambiguity of NNLL order, which corresponds to 
the intrinsic error of the NLL approximation. 

We present the results in the factorized form (15.180 . using the notation 



(D.l) 



Mx 



NLL 



Mo,xiQ') fr fx fl 



(D.2) 
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where X denotes a specific process. In particular we separate the finite parts of the correc- 
tions, and fx, from the subtracted electromagnetic part f^. The latter contains all 
soft /coUinear 1/e poles that must be cancelled against real photon emission or, in the case 
of initial-state singularities, factorized. In and fx we will omit contributions of 0{e) 
and C(e^), since such terms are irrelevant after cancellation of the photonic soft/collinear 
singularities. The electromagnetic contributions f^ have the general form 



a. 



^em NLL ^ ^ ,em ^ ( 

■'^ 4:71 ■'^'^ \4n 



2 rl , .2 



em 
2,X 



(D.3) 



with 



^^em NLL 



2e" 



L 



2 2 



2 6 



Mr 



2{L-l-k)Cl 
1 



cm 



^ad,em 
^l,X 



3 



+ (^6-3 + 2Le-2 + L^e-i) fegi^o} C^^lo + (2^^-^ + AL' - 

Crx,t + 0{e), (D.4) 



Le~^ + 2L^e-^ + 2L 



3\ 

"qed 



and the factors Cf^ ^, Cf'^ ^, Cf*^^™, which depend on the charges and masses of the 
external particles, are given in the next sections. The values for the /3-function coefficients 
6^^-* and 6qed can be found in (14. 5 p and (14. 6p . 



D.l Neutral-current four-fermion scattering 

This section deals with s-channel neutral-current four-fermion amplitudes 

Fk _pK^ ^k' 

J (J J(T J a' J a' 1 



(D.5) 



where a fermion-antifermion pair annihilates and produces another fermion-antifermion 
pair. The Born diagram of such an amplitude is given by 




We allow for top quarks only in the final state. Both particles of the initial state must 
share the same chirality k, and both particles of the final state must have the same 
chirality k', otherwise the amplitude is suppressed in NLL accuracy. The electromagnetic 
charge quantum numbers of the external particles are given hj qf = qf^ = —qj^ and 
g/' = = -gj.;, the hypercharges hy yj = yj. = -yj. and yf> = y^^' = 

and the isospin by 



the isospin components by tj 



and tj/ 



/3 



-/3 
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tf = tf/ = |. We use zj and zj, to denote the Yukawa factors (14.111) of the initial- 
and final-state particles, respectively. 

The amplitude can be written in the factorized form flS.lSp . 



A^NC — -^o,Nc(*5^) /nc /nc /nc • 



(D.6) 



The Born amplitude combined with the (non-diagonal) symmetric-electroweak contribu- 
tion /^"(^^ (Km reads 



47r 



47r/ 



Co,NC ^(2"^^ ~ ^"^^) (^i!nc) + L^^t ^ (2:/ + zj'j C'^NC + 5'2!Ncj 

+ 0(e), (D.7) 



_ r 3 /^ad /^sew 
'-^ '-'1,NC '-'l.NC 



where the Born term 



(D.8) 



is written in terms of couplings giiQ"^) renormalized at the scale Q, whereas and the 
other couplings and mixing angles in the loop corrections are renormalized at the scale /xr. 
The remaining terms in (ID. 71) read 



-2 „2 



4 



CTnc = %^ + §^/(^/ + !) + (/-/') 



ad 
NC 



Co, 



NC 



41n(^U^^ + 444'l-21nfe)Cf 



4 



sew 
2 / "1,NC 



sew 
y2,NC 



NLL /I ^3 , j2 

— Vr-'^ 



U 



, (1) ^7? V) , , (1) ^72 



.3.3 



+ 0{e). (D.9) 



The values for the /3-function coefficients 62^'' are given in (14.51) . The symmetric- 
electroweak result flD.7p is multiplied with the diagonal factors 

2 



rZ NLL . , , 

/nc - 4^ ^ 



^wt}-^„?) +(/"/') 



2 y 



0(e) 



(D.IO) 



and The latter is obtained from (lD.3p - (lD.4p with 



CI 



NC,0 



NC 



41n(-Jg/g/' -21n 2 



(D.ll) 
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where qt = 2/3, and the symbols 



1, = t 
0, ^ t 



f'fi 



/',t 



(D.12) 



are used to distinguish between massive and massless fermions in the final state. 

Note that only the electromagnetic contributions depend on the fermion masses. 
The remaining contributions are simply complemented by the Yukawa contributions pro- 
portional to in (lD.7p . otherwise they are equal to our results for massless fermions 
presented in Sect. 8.4.1 of Ref. [32]. 

D.2 Charged-current four-fermion scattering 

In this section we treat s-channel charged-current four-fermion amplitudes — > 
fa' fpi 1 where the fermions and f^' are the isospin partners of fp and fpi , respectively. 
We allow for a top-antibottom or bottom-antitop pair both in the initial and final state 
because the s-channel amplitude t b — t b arises via crossing symmetry as a contribution 
to the process b b ^ 1 1. 

For s-channel amplitudes with purely left-handed external fermions, 



J a J p J a' J p' 1 

the Born diagram involves the exchange of a W boson: 

fa\ / fa' 



(D.13) 




Most of the other combinations of chiralities for the external fermions yield contributions 
which are suppressed in NLL accuracy. The only exception is the s-channel amplitude 



t 



t^b^ 



(D.14) 



or, via crossing symmetry, b t — > b t , where the exchange of a (p scalar boson in the 
Born diagram produces a non-suppressed contribution with right-handed top quarks: 





We start with the fully left-handed amplitude (1D.13|) . The hypercharge quantum 
numbers of the external particles are given by yj = yjL = —yfL and ?/// = = — y/L, 

the isospin components by 



J a 



/3 

t- fL 
Jp 



/3 



and for left-handed fermions 
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= 1/2. The Yukawa factors (14. lip of the initial- and final-state particles are denoted 
by zj = zTl = zX and zj, = z\ = z\ , respectively. 

As in the previous section, the amplitude (1D.13|) is written in the form 



J\4qq — .Mo cciQ^) fee fee fee ■ 



2\ rsew cZ cem 



(D.15) 



Combining the Born amplitude with the (non-diagonal) symmetric-electroweak contribu- 
tion /^-(^ (Km . we find 



Mo,ceiQ') fee - Hp2, ^Yuip^, L) ui-p,, L^A'P^^ L) 



+ 



47r 



9im 



- 3L) Cr:Ee + L,^{zJ + zJ,)]+L 



"(ad 
,CC 



)2 



r 3 /~ia,d /^scw 



0{t 



(D.16) 



with 



CI 



sew 

cc 



'-'l.CC 



2 e2 ' 

41n(^^4M£:_2flnr^Ulnf^)14-21nf-^ 



tj 



CI 



sew NLL 



The Mz-dependent correction factor reads 



f'ee "i^^ 1 + ^ L f #4 + 2 ^4 ^) + 0[e) 



47r -^^^e^ ™ e2 " 4 
and the electromagnetic factor /q™ is obtained from (]D.3p - (ID.4p with 



72 4 + 7 72 + o-^/.t ?b + (/ ^ / ) 



4 e 



/^ad,em 
•-^1,00 



In 



-u 
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sew 
CC 



(D.17) 



(D.18) 



2 w 



tj e 



o2 w 



/ \ e2 * 



2 e2^/ ' 



(D.19) 



where gt = 2/3, gb = —1/3, and the symbols 



1, = t or fp = t 
0, otherwise 



1-5 



(D.20) 



and similarly for / — » /', /o- /o-/, /p — > /p/ are used to distinguish between massive and 
massless fermions. 

Again only the dependence of the electromagnetic contributions /^c on the fermion 
masses and the Yukawa terms in flD.16P are new compared with our results for massless 
fermions in Sect. 8.4.2 of Ref. [32]. 

Now we present the results for the s-channel amplitude t^b^ —>■ t^h^ f lD.14p which 
contributes via crossing symmetry to the process b'" b'" t^ t^. We need the hypercharges 
Utn = 4/3, ?/bL = 1/3 and the Yukawa factors 2;^ = 2, = 1. The amplitude reads 



•^CCY — -^CCCyIQ^) /coy /cCY /cCY 

and is expressed through 



(D.21) 



-^0,CCy(Q^) fcCY — 
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+ '4^ 
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ad 
CCY 



A? 



,CCY J + ^^^t ^ (^tR + ^bLj C^T^CCY + 92%CY 



T 3 y^ad y^scw 
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(D.22) 



with 



CI 



sew 
CCY 



CI 



ad 



CCY 



e2 4 4 e2 ' 

41n(^^^^-21nf^')cr 



,CCY 5 



NLL 



5'2,CCY I 2 

and the Mz-dependent factor 



ii' - ((,;■' I ^^i^^ + \ f) + (D.23) 
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The electromagnetic correction factor /^cy obtained from (]D.3p - (]D.4p with 

'-'1,CCY,0 ~ Vb' '-^l,' 



"Yem 
l,CCY,t 



CI 



ad, em 
CCY 



41n(-'lgtgb-21nf^) (ql + 



(D.25) 



D.3 Annihilation of two gluons into a fermion pair 

This section completes the four-particle processes by amplitudes with two gluons in 
the initial state. The process 



involves three diagrams at Born level, 



+ 




+ 



and the corresponding Born amplitude in the high-energy limit reads 



M 



-9s ^t^{Viyv{P2) U{-P3, k') X 

[ r'^t' [g'-'iPi - + g^'iPi + 2p2r - g'^c^Pi + P2y 



7p 



(D.26) 




u 



(D.27) 



where is the strong coupling, £^(^1,2) are the polarization vectors of the initial-state 
gluons with colour indices a and 6, are the generators of the QCD SU(3) gauge group, 
and f"'^'^ are the corresponding structure constants. 

As discussed in Sect. |6l our NLL results for n-fermion processes can be applied also to 
QCD processes involving fermions and gluons. The NLL amplitude for the process flD.261) 
assumes the usual factorized form. 



(D.28) 



and the presence of the gluons affects only the factorized Born amplitude A^o,g/- The 
NLL correction factors fgf", fgj, and /|™ are obtained from the general results of Sect. H] 
and Sect. Oby treating the reaction (ID. 260 as a — >■ 2 process, in the sense that the NLL 
correction factors receive contributions only from the final-state fermions. 
The symmetric-electroweak operator f^y is diagonal for this process. 



^sew NLL 



+ 



47r 
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,sew 
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(D.29) 
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and simply multiplies the Bom amplitude (1D.27|) . The terms Cf^^j- and gl™j- used here as 
well as the two other correction factors f^j and /|™ can be obtained from the corresponding 
results (ID.9p - (]D.lip of the neutral-current four-fermion amphtude in Sect. ID.ll by setting 
all electroweak quantum numbers of the initial-state fermions (?//, t^, tf, qj) to zero, 

If,, 



keeping only the final-state quantum numbers ?//', t?/, t//, 



D.4 Comparison with effective field-theory results 

The four-fermion amplitudes presented in this appendix can be compared with the 
results of Refs. [24,25] based on soft-collinear effective theory (SCET). To this end we 
have to use the results in Sect. VII of Ref. [25], omitting QCD contributions. We have 
found agreement at NLL accuracy for the symmetric-electroweak and Mz-dependent parts 
of our results. 

In the SCET framework, the full electroweak theory is matched at the high scale fi = Q 
to an effective theory SCETew where the degrees of freedom above the scale Q are inte- 
grated out. In NLL accuracy, only the tree-level expressions of the corresponding matching 
coefficients are relevant. These matching coefficients are evolved from the scale n = Q 
down to the scale /i = Mw using anomalous-dimension matrices calculated in SCETew- 
This step yields LL and NLL contributions which agree with the symmetric-electroweak 
parts /nc (1D-7P and f^Q (1D.16P of the neutral- and charged- current amplitudes presented 
here. At the low scale /x = Mw, SCETew is matched to another effective theory SCET^ 
where the massive gauge bosons are integrated out. The loop corrections resulting from 
this second matching agree with the factors /^c (ID.lOP and Jqq. flD.lSP arising in our 
calculation from the difference in the W- and Z-boson masses. 

Finally, the matching coefficients in SCET^ are evolved down to some finite scale 
yUo < Mw (/io = 30 GeV for the numerics in Ref. [25]), which acts as a cut-off for the 
singular contributions due to soft and coUinear photons. Due to the different regulariza- 
tion employed in our calculation these contributions cannot directly be compared with 
our electromagnetic factors /^c and /q™. 

In addition to the results of Ref. [25], our analysis in Sect. ID.21 also includes the 
charged- current t-channel amplitude for the process bb — > tt and, in particular, the 
tree-level exchange of a scalar boson in the case of right-handed top quarks. 

Our results concerning four-particle processes with two gluons and two quarks in 
Sect. ID. 31 are in agreement with the comments on these reactions in Ref. [25]. 
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